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Abstract 
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1. INTRODUCTION 



1. INTRODUCTION 

The code HFODD (v2.08k) solves the nuclear Skyrme-Hartree-Fock or Skyrme-Hartree-Fock- 
Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. It is a result of 
a long-term project that is developed at the Institute of Theoretical Physics, Warsaw University 
and Institut de Recherches Subatomiques, Universite Louis Pasteur in Strasbourg. The code 
has been published and is available from the Computer Physics Communications Program 



Library (http://www.cpc.cs.qub.ac.uk/cpc/ ). The description of the code and how to use it 
have been published in the Computer Physics Communications in the form of five independent 
publications, which are below referred to as I P, II |2j. Ill IV and V ISj. 

The present user's guide summarizes articles I-IV and provides practical details presented 
in the form that is suitable as a rapid reference on how to run the code. Information pertaining 
to the physics issues and numerical methods is not repeated here, and the user is referred to 
the original publications and to the references given therein. 

Together with the publication of the present user's guide we have established the HFODD 
home page (http:/ /www.fuw. edu.pl/~dobaczew/hfodd/hfodd. html), where more information is 
available, and which will contain information on future versions and extensions being developed. 



2. INPUT DATA FILE 
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2. INPUT DATA FILE 

The code reads the input data from the standard FORTRAN input file, and a general structure 
of the input data file is defined by the following rules: 

1. Input data file is an ASCII file composed of independent input items. 

2. An item is composed of the keyword line, and of the data line which follows immediately 
the keyword line. Two items contain no data lines, see point IHl below. 

3. The keyword line contains in columns 1 through 10 a keyword which is a specific text 
defining the item. If the keyword has less than 10 characters it has to be obligatorily 
padded with appropriate number of trailing spaces. Columns beyond 10 are ignored and 
can be used to place arbitrary comments or texts. 

4. The data line contains a prescribed number of data values which are read in the free 
FORTRAN FORMAT. The type of data (REAL, INTEGER, or CHARACTER) should match 
the definition of the given data line. The CHARACTER data should be provided starting 
from the 13-th column of the data line. 

5. Items can be separated by an arbitrary numbers of comment lines. A line is treated as 
a comment and ignored provided it does not contain in columns 1 through 10 any of the 
valid keywords. 

6. Two items contain no data lines. The first one is defined by the keyword EXECUTE, and 
starts a calculation for the currently defined input parameters. The second one is defined 
by the keyword ALL_DONE, and terminates the program. The latter item is not required 
if the compiler is able to properly recognize the end of input data file. When the end of 
input data file is encountered, the program terminates as if the keyword ALL_DONE was 
found. 

7. A given item may occur any number of times in different places of the input data file. 
Only the last one of the same items occurring before a given EXECUTE keyword is taken 
into account. 

8. The items appearing between two consecutive items EXECUTE define the data set, i.e, the 
set of input parameters, for which the calculation will start at the moment when the 
latter of the two such items is found. Between two consecutive items EXECUTE, the order 
of other items is arbitrary. 

9. All input parameters have the default values defined in the code (subroutine predef). 
Therefore, if the input data file contains only the item EXECUTE, and no other lines, a 
calculation for the default values is performed. In the present version, this results in 
performing 50 iterations for the superdeformed state of ^^^Dy at the angular frequency of 
hiu=0.5MeV. 

10. Only those items which define the values of input parameters which are different from the 
default values have to be included in the input data file. The values of input parameters 
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2. INPUT DATA FILE 



defined by any item stay in effect till another occurrence of the same item, or till the code 
terminates. This rule is valid irrespective of how many EXECUTE items follow the given 
item. 

In the following subsections we list all the valid keywords and their corresponding data lines. 
Titles of the subsections are introduced only to organize the keywords into groups pertaining 
to specific categories, but are otherwise irrelevant. The default values are given in data lines 
following the keyword lines. In the data lines we also show names of the corresponding variables, 
which are referred to in the text. 



2.1. General data 
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2.1. General data 

This section lists keywords pertaining to definitions of particle numbers, iterations, and phase 
space. 

2.1..1 Keyword: NUCLIDE 

86, 66 = IN_FIX, IZ_FIX 

Numbers of neutrons (IN_FIX) and protons (IZ_FIX) in the nucleus under consideration. Cal- 
culations for odd and odd-odd nuclei can be done only either with IR0TAT=1 or with IPAIRI=1 
and IPAHFB=0. In the latter case, the odd and odd-odd nuclei are described within the BCS 
blocking approximation with the time-odd potentials neglected. 



2.1.2. Keyword: ITERATIONS 

50 = NO ITER 

For NOITER>0, the specified number of iterations is performed. Unless IC0NTI=1 or IREAWS=1, 
the code starts with iteration number by constructing the initial mean fields based on the 
Nilsson potential. Specifying NOITER=0 requests only this initial phase. If the iteration is 
restarted from the previously recorded potentials (IC0NTI=1), then the counting of iterations 
continues from the number of iterations performed in the previous run. 



2.1.3. Keyword: ITERAT_EPS 

0.0 = EPSITE 



2.1.4. Keyword: MAXANTIOSC 

1 = NULAST 

These two parameters govern the termination of the HF iterations according to the achieved 
stability of solutions. The stability of the HF energy has been defined in (1-37) as the difference 
between the total energies calculated from the single-particle energies and from the Skyrme 
functional. The HF iterations continue until the absolute value of the stability is smaller than 
EPSITE (in MeV) over NULAST consecutive iterations. When this condition is fulfilled, iteration 
procedure terminates and the final results are printed. This allows for an automated adjustment 
of the number of iterations that are required to achieve a given level of convergence. The number 
of iterations NOITER, can now be set to a large value at which the iterations terminate if a stable 
solution is not found. 

The default value of EPSITE=0.0 deactivates this option, and then the code continues up 
to NOITER iterations regardless of the achieved stability. If a non-zero value of EPSITE is used, 
a non-zero value of NULAST has to be used too. In practice, a value of NULAST=5 prevents 
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2. INPUT DATA FILE 



an accidental termination of iterations in all cases when the stability energy gradually goes 
through zero, but the solution is not yet self-consistent. 



2.1.5. Keyword: PING-PONG 



0.0, 3 = EPSPNG, NUPING 



The code is able to detect the "ping-pong" divergence described in Sec. III-2.6, i.e., the situation 
when the HF iteration procedure gives oscillating results in every second iteration. Upon 
continuing the iteration, both sequences of results, i.e., those which correspond to the iteration 
numbers being even and odd, stay different but perfectly stable, and hence the correct self- 
consistent solution is never attained. 

The code recognizes such a situation by calculating the averages and variances of the stability 
energy (1-37), separately in the even and in the odd sequences of results, over the last NUPING 
pairs of iterations. The "ping-pong" divergence condition occurs when both variances become 
a factor EPSPNG smaller than the absolute value of the difference of the corresponding averages, 
i.e., when 



A{SS) 



even < EPSPNG X \6£cvcn - S£odd\ 
< EPSPNG X lIEcvcn - S£odd\ 



where 



(NUPING-1 
S£n-2p I /NUPING, 
p=0 

(NUPING-1 
E ^^n-2p-l I /NUPING, 
p=0 



(la) 
(lb) 



(2a) 
(2b) 



and 



A(5^)c 



(NUPING-1 
c 
p=0 



1/2 



/NUPING, 

1/2 



(NUPING-1 ^ 
E (^^n-2p-l -5?odd)' I /NUPING. 



(3a) 
(3b) 



Here, n denotes the number of the last accomplished HF iteration, and SSn' denotes the stability 
energy (1-37) obtained in the ?7,'-th iteration. 

The default value of EPSPNG=0.0 deactivates this option, and then the code continues up 
to NO ITER iterations regardless of the "ping-pong" divergence. If a non-zero value of EPSPNG 
is used, a value of NUPING>1 has to be used too. In practice, values of EPSPNG=0.01 and 
NUPING=3 allow for an efficient detection of the "ping-pong" divergence condition. 

Upon discovering the "ping-pong" divergence, the HF iterations are terminated and a table 
of absolute values of maximum differences of single-particle observables between the two se- 
quences of iterations is printed, see Sec. III-2.6. These maximum differences are determined for 



2.1. General data 



11 



states in each of the charge-parity-signature, charge-simplex, parity, or charge blocks, depend- 
ing on the conserved symmetries, and separately for particle and hole states. Whenever such 
a maximum difference is found for a particle state and for a hole state with adjacent indices, 
such a pair is proposed as a candidate for the diabatic blocking calculation. 



2.1.6. Keyword: CHAOTIC 

= NUCHAO 

The code is able to detect the "chaotic" divergence which occurs when the HF iterations give 
results which chaotically vary from one iteration to another one. The code recognizes such 
a divergence by finding the local maxima M^, k=l,2,. . . , in the sequence of absolute values 
of the stability energies (1-37), obtained in the entire series of the HF iterations performed. 
The "chaotic" divergence condition occurs when the code finds NUCHAO positive differences 
Mfc— Mfc_i. When this condition occurs, iteration procedure terminates and the final results 
are printed. 

For NUCHAO=0 (the default value) the code does not check whether the "chaotic" divergence 
occurs or not. In practice, a value of NUCHA0=5 allows for an efficient detection of the "chaotic" 
divergence condition. However, for a small value of NUCHAO and a small value of EPSPNG, the 
"ping-pong" divergence can sometimes be mistaken for the "chaotic" divergence. If one is inter- 
ested in the diabatic-blocking data, printed after the "ping-pong" divergence, the recommended 
value of NUCHA0=5 should be increased to 10 or more. 



2.1.7. Keyword: PHASESPACE 

0, 0, 0, =NUMBSP(0,0), NUMBSP(1,0), 
NUMBSP(0,1), NUMBSP(1,1) 

Numbers of the lowest mean-field eigenstates which are kept after the diagonalization of the 
mean-field Hamiltonians in the four charge-simplex blocks: {s,q) = {+i,n), {—i,n), {+i,p), 
{—i,p). All other eigenstates are discarded. If any of these numbers is equal to zero (the 
default value), the code sets it equal to the number of neutrons IN_FIX (for q=n) or protons 
IZ_FIX (for q=p). When no symmetry is conserved, the numbers of states kept in charge blocks 
are determined by sums NUMBSP(0,0) + NUMBSP(1,0) and NUMBSP(0,1) + NUMBSP(1,1) 

For calculations without pairing, i.e. for IPAIRI=0, the user is responsible for using values of 
NUMBSP large enough to accommodate all wave functions which might be useful for the required 
vacuum and particle-hole configurations, see Sec, II-3.4. In practice, the use of the default 
values described above is recommended as a safe option. The size of the matrices defined by 
the NDSTAT parameter can be reduced by the user when small values of NUMBSP are used. 

For calculations with the BCS pairing, i.e. for IPAIRI=1 and IPAHFB=0, values of NUMBSP 
define the pairing window in which the BCS equations are solved. 

For calculations with the HFB pairing, i.e. for IPAIRI=1 and IPAHFB=1, values of NUMBSP 
define the numbers of canonical states kept after diagonalizing density matrices. They are 
irrelevant for the pairing window, which is then defined by the cutoff energy ECUTOF. 
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2.2. Interaction 

This section lists keywords pertaining to definitions of coupling constants and interaction pa- 
rameters. 

2.2.1. Keyword: SKYRME-SET 



CHARACTER*4 acronym of the Skyrme force parameter set. Must start at the 13-th column of 
the data line. In version (v2.08k), vahd acronyms are SIII, SV, SKM*, SKP, SKMP, SKIl, 
SKO, SKOP, SLY4, SLY5, MSKl, MSK2, MSK3, MSK4, MSK5, MSK6, SKX, SKXC. Other 
sets of parameters can easily be included in the subroutine params. 



2.2.2. Keyword: SKYRME-STD 



Parameters of several standard Skyrme forces are encoded within the program. Calculation for 
a given standard force can be requested by specifying in the input data file its acronym SKYRME. 
Along with the force parameters, for each force there is encoded information on how the given 
force should be used, i.e., with which value of the parameter ^^/2m, and with which treatment 
of tensor, spin-orbit, and center-of-mass terms. For ISTAND=1, calculations are performed with 
these features set in the way specific for the given force, and the rest of the switches read on the 
same line is ignored. For ISTAND=0, the switches define nonstandard features in the following 
way: 

• For KETA_J=1 the code takes into account in the functional, and for KETA_J=0 neglects, 
the so-called tensor terms. The second option is equivalent to setting in Eq. (1-12) 
coupling constants C/=0 and Cj=Q. 

• For KETA_W=0 or 1 the code uses the traditional {Wq = Wq) or generalized (Wq ^ Wq) 
spin-orbit term, respectively, where 



For forces that use generalized spin-orbit term, option KETA_W=0 sets equal to Wq^ 
while for forces that use traditional spin-orbit, option KETA_W=1 has no effect. For 
KETA_W=2, strengths Wq and are set equal to the values of WO_INP and WOPINP, re- 
spectively, which are read under the keyword of SPIN_ORBIT. Note that if KETA_W=2 is 
used, and keyword SPIN_ORBIT is not specified in the input data file, the default values 
of WO_INP and WOPINP supersede those that are encoded within the program for the given 
Skyrme force. 



SKM* = SKYRME 



0, 1, 0, 0, = ISTAND,KETA_J,KETA_W,KETACM,KETA_M 




(4a) 
(4b) 



2.2. Interaction 
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For KETACM=0 the code uses the traditional one-body center-of-mass correction before 
variation, 

Ee.m.^i?f;. = -^(T), (5) 

uses the two-body center-of-mass correction after variation 



where T=—^ Hi=i is the one-body kinetic-energy operator. For KETACM=1 the code 



where -P=— «^ I^iLi is the total hnear momentum operator. Value of KETACM=2 is 
reserved for a future implementation of the two-body center-of-mass correction before 
variation. For KETACM=3 the center-of-mass correction is neglected. 

For KETA_M=0 the code uses the value of ^^/2m=20. 73620941 MeVfm^, which was encoded 
in version (vl.75r). For KETA_M=1 the code uses the value specific for the given Skyrme 
force. For KETA_M=2 the code uses the value specified in the input data file under keyword 
HBAR20VR2M. 



2.2.3. Keyword: HBAR20VR2M 

20.73620941 = HBMINP 

Value of the Ti^ /2m parameter, which is used if KETA_M=2 is set under keyword SKYRME-STD, 
and ignored otherwise. 



2.2.4. Keyword: SPIN_ORBIT 

120.0, 120.0 = WO_INP,WOPINP 

Strengths Wq and Wq of the generalized spin-orbit interaction (@]), which are used if KETA_W=2 
is set under keyword SKYRME-STD, and ignored otherwise. 



2.2.5. Keyword: LANDAU 

0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 =LAN0DD,X0_LAN,X1_LAN, 

GO_LAN,GOPLAN, 
G1_LAN,G1PLAN 

Three-digit steering switch LANODD is used to specify the type of operations performed in order 
to define selected time-odd coupling constants based on values of the Landau parameters. The 
switch is followed by values of xq, xi, g^, q'q, gi, and g[, as given by equations: 

xo = C^[0]/C^[ps.t], (7a) 
xi = Cm/Cnps.t], (7b) 
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N,{2C^,[p,,,]+2Cl(3plil), (7c) 
g', = Noi2C([p.,,]+2Cf/3pli^), (7d) 
g, = -2N,C^(3plLl (7e) 



g[ = -2N,C^PpILI (7f) 




where {3 = (37rV2)2/3 and 



For LANODD=000, coupling constants are not defined by the Landau parameters, and the rest 
of parameters read on the same data hne is ignored. For LAN0DD=111, Eqs. (|7j) are solved, and 
coupling constants Co[0], Co[psat]) C*i[0], C'liPsat], Cq, and Cf are determined from the values 
of Xo, Xi, go, (7q, gi, and g[. For other values of LANODD, the following three steps are performed 
in sequence: 

1. Step one: If the rightmost digit of LANODD is equal to 1 then the coupling constants Cq 
and Cf are determined form gi and g[, Eqs. (|7e|) and (ffH): otherwise these two coupling 
constants are determined from the Skyrme force parameters. 

2. Step two: If the middle digit of LANODD is equal to 1 then the coupling constants CQ[psat] 
and Cf[psat] are determined form go and g'o, Eqs. (fTcjl and (|7d|l : otherwise these two 
coupling constants are determined from the Skyrme force parameters. 

3. Step three: If the leftmost digit of LANODD is equal to 1 then the coupling constants Co[0] 
and C([0] are determined form xq and xi, Eqs. (f7a|) and (f7b|) : otherwise these two coupling 
constants are determined from the Skyrme force parameters. 

All combinations of zeros and ones are allowed in LANODD. 



2.2.6. Keyword: LANDAU-SAT 

-1.0, -LO, -1.0 = HBMSAT,RHOSAT,EFFSAT 

Values of parameters /2m, m* /m, and psat, respectively, which are used when solving Eq. 
((Tj). If a negative number is read for any of these parameters, then the program uses the 
corresponding value calculated from parameters of the given Skyrme force. 



2.2.7. Keyword: EVE_SCA_TS 



1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 

SRHO_T, SRHO_S,SRHODT, SRHODS,SLPR_T, SLPR_S, 

STAU_T, STAU_S, 
SSCU_T, SSCU_S, 
SDIV_T, SDIV_S 



By using this item, the coupling constants corresponding to a given Skyrme parameter set can 
be arbitrarily scaled. This allows calculations with modified Skyrme functionals. The time-even 
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coupling constants in the total-sum representation, 



'-"tot 



r<P r^P 
•-"O "~ '-"1; 

2Cf, 



(9a) 
(9b) 



are multiplied by the 12 numbers from SRHO_T to SDIV_S. Variables with names ending with _T 
and _S multiply the "total" and "sum" coupling constants, respectively. Variables with names 
containing the acronyms RHO, LPR, TAU, SCU, and DIV, multiply the coupling constants with 
superscripts p, Ap, r, J, and VJ, respectively, and those with RHOD multiply the density- 
dependent part of C . Similar name convention is used for many other variables in the code 

HFODD. 



2.2.8. Keyword: ODD_SCA_TS 



Same as above but for the time-odd coupling constants. Acronyms SPI, LPS, CUR, KIS, and 
ROT correspond to coupling constants with superscripts s. As, T, j, and Vj, respectively, and 
those with SPID correspond to the density-dependent part of C"^. 



2.2.9. Keyword: EVE_SCA_PM 



Same as above but for the time-even coupling constants in the isoscalar-isovector representation. 



1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 

SSPI_T, SSPI_S,SSPIDT, SSPIDS,SLPS_T, SLPS_S, 



SCUR_T, SCUR_S, 
SKIS_T, SKIS_S, 
SROT_T, SROT_S 



1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 

SRHO_P, SRHO_M,SRHODP, SRHODM,SLPR_P, SLPR_M, 



STAU_P, STAU_M, 
SSCU_P, SSCU_M, 
SDIV_P, SDIV_M 




tot ) 



(10a) 
(10b) 



Variables with names ending with _P and _M multiply the isoscalar and isovector coupling 
constants, respectively. The total-sum scaling factors are used first, and the isoscalar-isovector 
scaling factors are used afterwards. 
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2. INPUT DATA FILE 



2.2.10. Keyword: ODD_SCA_PM 

1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 

SSPI_P, SSPI_M,SSPIDP, SSPIDM,SLPS_P, SLPS_M, 

SCUR_P, SCUR_M, 
SKIS_P, SKIS_M, 
SROT_P, SROT_M 

Same as above but for the time-odd coupling constants in the isoscalar-isovector representation. 



2.2.11. Keyword: G_SCALING 

1.0, 1.0 = FACTGN, FACTGP 

For IPAIRI=1 the code HFODD solves the BCS equations with the neutron and proton pairing 
strengths defined in Ref. [H]. These values can be modified by defining here the multiplicative 
factors FACTGN and FACTGP for neutrons and protons, respectively. 



2.2.12. Keyword: INI_FERMI 

-8.0, -8.0 = FERINI(0),FERINI(1) 



2.2.13. Keyword: INI_DELTA 

1.0, 1.0 = DELINI(0),DELINI(1) 

For IPCONT=0, the HFB calculations that are restarted from previously saved results are per- 
formed with new values of the neutron and proton Fermi energies, FERINI(O) and FERINI(l), 
and pairing gaps, DELINI(O) and DELINK 1). New values of pairing gaps are implemented by 
overwriting the old HFB pairing potentials with constant values of DELINI(O) and DELINK 1) 
for neutrons and protons, respectively. These constant potentials are ignored after the first 
iteration, i.e., in the first iteration, mixing of previous and current potentials (see keyword 
SLOW-PAIR) is not performed. A possibility of restarting calculations with nonzero pairing is 
very useful in case the pairing would have vanished in a former run. For IPAHFB=0 or IPCONT=0, 
values of FERINI and DELINI are ignored. 



2.2.14. Keyword: FIXDELTA_M 

1.0, = DELFIN,IDEFIN 

For IDEFIN=1, the HFB pairing calculations are performed with a fixed value of the neutron 
pairing gap equal to DELFIN. For IDEFIN=0, value of DELFIN is ignored and the zero-range 
density-dependent pairing force (fTT|) is used for neutrons. IDEFIN=1 requires IPAHFB=1. 
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2.2.15. Keyword: FIXDELTA_P 

1.0, = DELFIP,IDEFIP 

Same as above but for the proton pairing gap. 



2.2.16. Keyword: PAIRNFORCE 

-200.0, 0.0, 1.0 = PRHO_N,PRHODN,POWERN 

Parameters Vq, Vi, and a, respectively, of the zero-range density-dependent pairing force defined 
by the form-factor: 

fir) = Vo + yip"(r) =Vo(^l- (^^) ) , (11) 

which is used in the HFB calculations for neutrons. In case values of Vi and a allow it, the 
code calculates the value of the saturation density po that gives the equivalent form of the 
form-factor (fTT|) . In case po cannot be calculated, the codes set its value to 1; po is calculated 
only for the purpose of information, while internally the code uses only the value of Vi. For 
IDEFIN=1 or IPAHFB=0, parameters PRHO_N,PRHODN,POWERN are ignored. 



2.2.17. Keyword: PAIRPFORCE 

-200.0, 0.0, 1.0 = PRHO_P,PRHODP,POWERP 

Same as above but for the proton pairing force. 



2.2.18. Keyword: PAIR_FORCE 

-200.0, 0.0, 1.0 = PRHO_T,PRHODT,POWERT 

Same as above but for the neutron and proton pairing force. This keyword is equivalent to 
using two keywords simultaneously, PAIRNFORCE and PAIRPFORCE, with identical parameters 
for neutrons and protons. 



2.2.19. Keyword: PAIRNINTER 

-200.0, 0.16, 1.0 = PRHO_N,PRHOSN,POWERN 

Parameters Vq, po, and a, respectively, of the zero-range density-dependent pairing force (fTT|) 
used in the HFB calculations for neutrons. In case values of po and a allow it, the code calculates 
the value of Vi that gives the equivalent form of the form-factor ()11|) . In case Vi cannot be 
calculated, the code stops. For IDEFIN=1 or IPAHFB=0, parameters PRHO_N,PRHOSN,POWERN are 
ignored. 
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2. INPUT DATA FILE 



2.2.20. Keyword: PAIRPINTER 

-200.0, 0.16, 1.0 = PRHO_P,PRHOSP,POWERP 

Same as above but for the proton pairing force. 

2.2.21. Keyword: PAIR_INTER 

-200.0, 0.16, 1.0 = PRHO_T,PRHOST,POWERT 

Same as above but for the neutron and proton pairing force. This keyword is equivalent to 
using two keywords simultaneously, PAIRNINTER and PAIRPINTER, with identical parameters 
for neutrons and protons. 



2.2.22. Keyword: CUTOFF 

60.0 = ECUTOF 

The cutoff energy Cmax for summing up contributions of the HFB quasiparticle states to density 
matrices, see Sec. IV-2.5. Ignored for IPAHFB=0. 
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2.3. Symmetries 

This section lists keywords pertaining to selection of conserved and broken symmetries. 
2.3.1. Keyword: SIMPLEXY 

1 = ISIMPY 

For ISIMPY=1, calculation with y-simplex conserved are performed, while for ISIMPY=0 the 
simplex is broken. Value of ISIMPY must be consistent with switches IPARTY, ISIGNY, ISIMTY, 
and IROTAT, see Tables □ and H 



2.3.2. Keyword: SIGNATUREY 

1 = ISIGNY 

For ISIGNY=1, calculation with signature conserved are performed, while for ISIGNY=0 the 
signature is broken. Value of ISIGNY must be consistent with switches ISIMPY, IPARTY, ISIMTX, 
and ISIMTZ, see Tables □ and El 



2.3.3. Keyword: PARITY 

-1 = IPARTY 

For IPARTY=1, calculation with parity conserved are performed, while for IPARTY=0 the parity 
is broken. Value of IPARTY must be consistent with switches ISIMPY and ISIGNY, see Table 
lU For IPARTY=-1 (the compatibility mode), the code sets IPARTY=ISIGNY and requires that 
ISIMPY=1 



2.3.4. Keyword: ROTATION 

1 = IROTAT 

For IR0TAT=1, calculation with time-reversal breaking are performed, while for IROTAT=0 the 
time-reversal symmetry is conserved. In the latter case the calculations are performed only 
for one value of the simplex, s=+i, which gives almost twice shorter CPU times. IROTAT=0 
is incompatible with providing a non-zero value of the angular frequency or with attempting a 
calculation for an odd or odd-odd nucleus. IR0TAT=1 is incompatible with IPAIRI=1, unless 
IPAHFB=1, i.e., the rotating solutions can be obtained for the HFB pairing but not for the BCS 
pairing. 



2.3.5. Keyword: TIMEREVERS 

= ITIREV 
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2. INPUT DATA FILE 



Table 1: Primary set of conserved and nonconserved symmetries allowed in the code hfodd 
version (v2.08k). 



Option 


Symmetries 


Switches 


Sy 


Ry 


P 


ISIMPY 


ISIGNY 


IPARTY 


PI 


conserved 


conserved 


conserved 


1 


1 


1 


P2 


conserved 


nonconserved 


nonconserved 


1 








P3 


nonconserved 


conserved 


nonconserved 





1 





P4 


nonconserved 


nonconserved 


conserved 








1 


P5 


nonconserved 


nonconserved 


nonconserved 












Table 2: Secondary set of conserved and nonconserved symmetries allowed in the code hfodd 
version (v2.08k). 



Option 


Symmetries 


Switches 


Sy 


f 


^y 


ISIMPY 


ITIREV 


ISIMTY 


SI 


conserved 


conserved 


conserved 


1 


1 


1 


S2 


conserved 


nonconserved 


nonconserved 


1 








S3 


nonconserved 


conserved 


nonconserved 





1 





S4 


nonconserved 


nonconserved 


conserved 








1 


S5 


nonconserved 


nonconserved 


nonconserved 












Table 3: Ternary set of conserved and nonconserved symmetries allowed in the code hfodd 
version (v2.08k). 



Option 


Symmetries 


Switches 


Ry 


si 


qT 


ISIGNY 


ISIMTX 


ISIMTZ 


Tl 


conserved 


conserved 


conserved 


1 


1 


1 


T2 


conserved 


nonconserved 


nonconserved 


1 








T3 


nonconserved 


conserved 


nonconserved 





1 





T4 


nonconserved 


nonconserved 


conserved 








1 


T5 


nonconserved 


nonconserved 


nonconserved 
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For ITIREV=1, calculation with time-reversal conserved are performed, while for ITIREV=0 this 
symmetry is broken. This switch is used only as a convenient replacement of switch IROTAT; 
ITIREV=1 is equivalent to IROTAT=0 and ITIREV=0 is equivalent to IR0TAT=1. 



2.3.6. Keyword: TSIMPLEX_Y 

-1 = ISIMTY 

For ISIMTY=1, calculation with y-simplex"^ conserved are performed, while for ISIMTY=0 this 
symmetry is broken. Value of ISIMTY must be consistent with switches ISIMPY and IROTAT, 
see Table El For ISIMTY=-1 (the compatibility mode), the code sets ISIMTY to 1-IROTAT and 
requires that ISIMPY=1. 



2.3.7. Keyword: TSIMPLEXES 

1, 1 = ISIMTX, ISIMTZ 

For ISIMTX=1 and/or ISIMTZ=1, calculation with conserved symmetries given by x-simplex^ 
and/or z-simplex^ are performed, respectively. These symmetries are broken for ISIMTX=0 and 
ISIMTZ=0. Values of ISIMTX and ISIMTZ must be consistent with that of ISIGNY, see Table El 



2.3.8. Keyword: TSIMPLEX3D 

1, -1, 1 = ISIMTX,ISIMTY,ISIMTZ 

This keyword allows to simultaneously input all the three switches that define the three T- 
simplexes. It is equivalent to using keywords TSIMPLEXES and TSIMPLEX_Y together. 



2.3.9. Keyword: PAIRING 

= IPAIRI 

For IPAIRI=0, calculation without pairing correlations (pure HF) are performed, while for 
IPAIRI=1 the pairing correlations are included. In the latter case, for IPAHFB=0 or 1 pairing 
is included within the BCS or HFB method, respectively. IPAIRI=1 is incompatible with 
IR0TAT=1, unless IPAHFB=1, i.e., the rotating solutions can be obtained for the HFB pairing 
but not for the BCS pairing. 



2.3.10. Keyword: HFB 

= IPAHFB 

For IPAHFB=1, calculation with pairing correlations included by solving the HFB equation are 
performed. IPAHFB=1 requires IPAIRI=1. IPAHFB=1 also requires ISIMPY=1, i.e., for broken 
simplex symmetry, solution of the HFB equation is in version (v2.08k) not yet implemented. 
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2. INPUT DATA FILE 



2.4. Configurations with no conserved symmetries 

This section lists keywords pertaining to definitions of configurations in the case when no 
symmetries are conserved, and the single-particle or single-quasiparticle space is divided into 
two charge blocks (neutrons and protons). 



2.4.1. Keyword: PHNONE_NEU 

1, 00, 00 = NUPAHO,KPNONE,KHNONE 

Neutron particle-hole excitations pertaining to the situation when all neutrons are in one com- 
mon block (no simplex, signature, or parity is conserved). NUPAHO is the consecutive number 
from 1 to 5 (up to five sets of excitations can be specified in separate items). Particles are re- 
moved from the KHNONE-th state and put in the KPNONE-th state. At every stage of constructing 
excitations the Pauli exclusion principle has to be respected (particle removed from an occupied 
state and put in an empty state). Values equal zero have no effect. Note that for all neutrons 
sitting in one common block the reference configuration from which the particle-hole excita- 
tions are counted is defined by the total number of neutrons. These particle-hole excitations 
are ignored unless ISIMPY=0, IPARTY=0, and IPAIRI=0. 

2.4.2. Keyword: PHNONE_PRO 

1, 00, 00 = NUPAHO,KPNONE,KHNONE 
Same as above but for the proton particle-hole excitations. 

2.4.3. Keyword: DIANON_NEU 

2, 1, = KPFLIZ,KHFLIZ,KOFLIZ 

The diabatic blocking of neutron single-particle configurations in the situation when all neutrons 
are in one common block. The rules to define the blocking are analogous to those described for 
the conserved parity, keyword DIAPAR_NEU. 



2.4.4. Keyword: DIANON_PRO 

2, 1, = KPFLIZ,KHFLIZ,KOFLIZ 

Same as above but for the diabatic blocking of proton single-particle configurations in the 
situation when all protons are in one common block. 



2.5. Configurations with conserved simplex 
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2.5. Configurations with conserved simplex 

This section lists keywords pertaining to definitions of configurations in the case when simplex 
is conserved, and the single-particle or single-quasiparticle space is divided into four charge- 
simplex blocks (neutrons and protons with simplexes s=±i). 



2.5.1. Keyword: VACSIM_NEU 

43, 43 = KVASIM(O), KVASIM(l) 

Numbers of lowest neutron states occupied in the two simplex blocks, denoted by (+) and (— ), 
of given simplexes, s=+i and s=—i, respectively. These numbers define the simplex reference 
configuration from which the particle-hole excitations are counted. The definitions of simplex 
reference configuration and excitations are ignored unless ISIMPY=1, ISIGNY=0, and IPAIRI=0. 



2.5.2. Keyword: VACSIM_PRO 

33, 33 = KVASIM(O), KVASIM(l) 

Same as above but for the numbers of proton states. 

2.5.3. Keyword: PHSIMP_NEU 

1, 00, 00, 00, 00 = NUPAHO, KPSIMP, KPSIMM, KHSIMP, KHSIMM 

Neutron particle-hole excitations in the simplex blocks. NUPAHO is the consecutive number from 
1 to 5 (up to five sets of excitations can be specified in separate items). Particles are removed 
from the KHSIMP-th state in the (+) block and from the KHSIMM-th state in the (— ) block, 
and put in the KPSIMP-th state in the (-I-) block and in the KPSIMM-th state in the (— ) block. 
At every stage of constructing excitations the Pauli exclusion principle has to be respected 
(particle removed from an occupied state and put in an empty state). Values equal zero have 
no effect. In practice, reasonable excitations can only be constructed by consulting the printed 
lists of single-particle states with their consecutive numbers in given blocks. These particle-hole 
excitations are ignored unless ISIMPY=1, ISIGNY=0, and IPAIRI=0. 



2.5.4. Keyword: PHSIMP_PRO 

1, 00, 00, 00, 00 = NUPAHO, KPSIMP, KPSIMM, KHSIMP, KHSIMM 
Same as above but for the proton particle-hole excitations. 
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2.5.5. Keyword: DIASIM_NEU 

2,2, 1,1, 0,0= KPFLIM(0,0),KPFLIM(1,0), 

KHFLIM (0 , 0) ,KHFLIM (1,0), 
K0FLIM(0,0),K0FLIM(1,0) 

The diabatic blocking of neutron single-particle simplex configurations. Matrices KPFLIM con- 
tain the indices of the particle states in the two simplex blocks denoted by (+) and (— ), of given 
simplex values, i.e., s=+i and —i, respectively. Matrices KHFLIM contain analogous indices of 
the hole states, and matrices KOFLIM define types of blocking in analogy to those shown in 
Table H 



2.5.6. Keyword: DIASIM_PRO 

2,2, 1,1, 0,0= KPFLIM(0,1),KPFLIM(1,1), 

KHFLIM(0,1),KHFLIM(1,1), 
K0FLIM(0,1),K0FLIM(1,1) 



Same as above but for the diabatic blocking of proton single-particle simplex configurations. 



2.6. Configurations with conserved parity 



25 



2.6. Configurations with conserved parity 

This section lists keywords pertaining to definitions of configurations in the case when parity is 
conserved, and the single-particle or single-quasiparticle space is divided into four charge-parity 
blocks (neutrons and protons with parities 7r=±l). 



2.6.1. Keyword: VACPAR_NEU 

44, 42 = KVASIQ(0),KVASIQ(1) 

Numbers of lowest neutron states occupied in the two parity blocks, denoted by (+1) and 
(— 1), of given parities, 7r=+l and 7r=— 1, respectively. These numbers define the parity ref- 
erence configuration from which the particle-hole excitations are counted. The definitions of 
parity reference configuration and excitations are ignored unless IPARTY=1, or IPARTY=— 1 and 
ISIMPY=ISIGNY=1. They are also ignored for ISIMPY=1 and IPAIRI=1. 



2.6.2. Keyword: VACPAR_PRO 

32, 34 = KVASIQ(0),KVASIQ(1) 
Same as above but for the numbers of proton states. 



2.6.3. Keyword: PHPARI_NEU 

1, 00, 00, 00, 00 = NUPAHO,KPSIQP,KPSIQM,KHSIQP,KHSiqM 

Neutron particle-hole excitations in the parity blocks. NUPAHO is the consecutive number from 
1 to 5 (up to five sets of excitations can be specified in separate items). Particles are removed 
from the KHSIQP-th state in the {+1) block and from the KHSIQM-th state in the (—1) block, and 
put in the KPSIQP-th state in the (+1) block and in the KPSIQM-th state in the (-1) block. At 
every stage of constructing excitations the Pauli exclusion principle has to be respected (particle 
removed from an occupied state and put in an empty state). Values equal zero have no effect. 
These particle-hole excitations are ignored unless ISIMPY=0, IPARTY=1, and IPAIRI=0. 



2.6.4. Keyword: PHPARI_PRO 

1, 00, 00, 00, 00 = NUPAHO,KPSIQP,KPSIQM,KHSiqP,KHSIQM 

Same as above but for the proton particle-hole excitations. 
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2.6.5. Keyword: DIAPAR_NEU 

2,2, 1,1, 0,0= KPFLIQ(0,0),KPFLIQ(1,0), 

KHFLIQ (0 , 0) ,KHFLIQ (1 , 0) , 
K0FLIQ(0,0),K0FLIQ(1,0) 

The diabatic blocking of neutron single-particle parity configurations. Matrices KPFLIQ contain 
the indices of the particle states in the two parity blocks denoted by (+1) and (—1), of given 
parities, i.e., 7r=+l and —1, respectively. Matrices KHFLIQ contain analogous indices of the 
hole states, and matrices KOFLIQ define the type of blocking according to definitions shown in 
Table H 



Table 4: Values of parameter KOFLIQ that specify type of blocking in case of conserved parity. 
Analogous values define type of blocking for other conserved symmetries. 



KOFLiq=0 

KOFLiq=+l 

KOFLiq=-l 

K0FLiq=+2 

K0FLiq=-2 

K0FLiq=+3 

K0FLiq=-3 

K0FLiq=+4 

K0FLiq=-4 

K0FLIQ=+5 

K0FLiq=-5 

K0FLiq=+6 

K0FLiq=-6 

K0FLiq=+7 

K0FLiq=-7 

K0FLiq=+8 

K0FLiq=-8 

K0FLiq=+9 

K0FLiq=-9 

K0FLiq=+10 

K0FLiq=-10 



No diabatic blocking 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 
The state which has 



; in the given parity block. 

the larger ^/-alignment is occupied. 

the smaller ^/-alignment is occupied. 

the larger y-intrinsic spin is occupied. 

the smaller y-intrinsic spin is occupied. 

the larger x-alignment is occupied. 

the smaller x-alignment is occupied. 

the larger x-intrinsic spin is occupied. 

the smaller x-intrinsic spin is occupied. 

the larger z-alignment is occupied. 

the smaller ^-alignment is occupied. 

the larger z-intrinsic spin is occupied. 

the smaller ^-intrinsic spin is occupied. 

the larger o-alignment is occupied. 

the smaller o-alignment is occupied. 

the larger o-intrinsic spin is occupied. 

the smaller o-intrinsic spin is occupied. 

the larger multipole moment Q20 is occupied. 

the smaller multipole moment Q20 is occupied. 

the larger multipole moment Q22 is occupied. 

the smaller multipole moment Q22 is occupied. 



Here, the x-, y-, and ^-alignments or intrinsic spins denote projections of the total angular 
momentum or spin, respectively, on the x, y, and z axes. Similarly, the o-alignment or o-intrinsic 
spin denotes analogous projections on the direction of the angular frequency ujj. 

Within the diabatic blocking procedure one does not predefine whether the particle or the 
hole state is occupied (as is the case when the particle-hole excitations are defined, see Section 
3.4 of II). In each iteration the code calculates the average alignments (or average intrinsic 
spins, or average quadrupole moments) of both states (those defined by KPFLIM and KHFLIM), 
and occupies that state for which a larger, or a smaller value is obtained. Therefore, the order 



2.6. Configurations with conserved parity 
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of both states in the Routhian spectrum is irrelevant. 

The user is responsible for choosing the particle-state indices (in KPFLIM) only among those 
corresponding to empty single-particle states, and the hole-state indices (in KHFLIM) only among 
those corresponding to occupied single-particle states. 



2.6.6. Keyword: DIAPAR_PRO 

2,2, 1,1, 0,0= KPFLIQ(0,1),KPFLIQ(1,1), 

KHFLIQ (0 , 1) ,KHFLIQ (1 , 1) , 
KOFLIQ (0 , 1) ,KOFLIQ (1 , 1) 



Same as above but for the diabatic blocking of proton single-particle parity configurations. 
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2.7. Configurations with conserved parity and signature 

This section lists keywords pertaining to definitions of configurations in tlie case wfien parity 
and signature are conserved, and tlie single-particle or single-quasiparticle space is divided into 
eight charge-parity-signature blocks (neutrons and protons with parities 7r=±l and signatures 
r=±i). 



2.7.1. Keyword: VACSIG_NEU 

22, 22, 21, 21 =KVASIG(0,0), KVASIG(0,1), 
KVASIG(1,0), KVASIG(1,1) 

Numbers of lowest neutron states occupied in the four parity-signature blocks, denoted by 
(+,+), (+,—), (—,+), and (—,—), of given (parity,signature) combinations, i.e., (vr, r) = (+l,+i), 
(+1,-2), (— 1,+z), and (— 1,— i), respectively. These numbers define the parity-signature ref- 
erence configuration from which the particle-hole excitations are counted. The definitions of 
parity-signature reference configuration and excitations are ignored unless ISIMPY=1, ISIGNY= 
1, and IPAIRI=0. 



2.7.2. Keyword: VACSIG_PRO 



16, 16, 17, 17 =KVASIG(0,0), KVASIG(0,1), 
KVASIG(1,0), KVASIG(1,1) 



Same as above but for the numbers of proton states. 



2.7.3. Keyword: PHSIGN_NEU 

1, 00,00,00,00, 00,00,00,00 =NUPAHO, 

KPPPSP, KPPPSM, KPPMSP, KPPMSM, 
KHPPSP, KHPPSM, KHPMSP, KHPMSM 

Neutron particle-hole excitations in the parity-signature blocks. Basic principles are the same 
as for the excitations in the parity blocks. Particles are removed from the KHPPSP-th state in 
the (+,+) block, from the KHPPSM-th state in the (+,-) block, from the KHPMSP-th state in the 
(-,+) block, and from the KHPMSM-th state in the (-,-) block, and put in the KPPPSP-th state 
in the {+,+) block, in the KPPPSM-th state in the (+,-) block, in the KPPMSP-th state in the 
(— ,-|-) block, and in the KPPMSM-th state in the (— ,— ) block. These particle-hole excitations 
are ignored unless ISIMPY=1, ISIGNY=1, and IPAIRI=0. 
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2.7.4. Keyword: PHSIGN_PRO 

1, 00,00,00,00, 00,00,00,00 =NUPAHO, 

KPPPSP, KPPPSM, KPPMSP, KPPMSM, 
KHPPSP, KHPPSM, KHPMSP, KHPMSM 

Same as above but for the proton particle-hole excitations. 



2.7.5. Keyword: DIASIG_NEU 

2, 2, 2, 2, 1, 1, 1, 1, 0, 0, 0, = 

KPFLIG (0 ,0,0) ,KPFLIG (0 ,1,0) ,KPFLIG (1,0,0) ,KPFLIG ( 1,1,0), 
KHFLIG (0 ,0,0) ,KHFLIG (0 ,1,0) ,KHFLIG (1,0,0) ,KHFLIG ( 1,1,0), 
KOFLIG (0 ,0,0) ,KOFLIG (0 ,1,0) ,KOFLIG (1,0,0) ,KOFLIG ( 1,1,0) 

The diabatic blocking of neutron single-particle parity-signature configurations. Matrices 
KPFLIG contain the indices of particle states in the four parity-signature blocks denoted by 
(-|-,+), (+,—), (—,+), and (—,—), of given (paritysignature) combinations, i.e., {n,r) = (-1-1, +z), 
(-1-1,— z), (— 1,-fz), and respectively. Matrices KHFLIG contain analogous indices of hole 

states, and matrices KOFLIG define the type of blocking according to Table IH Other rules de- 
scribed for the parity case apply here analogously. 



2.7.6. Keyword: DIASIG_PRO 

2, 2, 2, 2, 1, 1, 1, 1, 0, 0, 0, = 

KPFLIG(0 , , 1) ,KPFLIG(0 ,1,1) ,KPFLIG (1,0, 1) ,KPFLIG( 1 ,1,1), 
KHFLIG (0 , , 1) , KHFLIG (0 ,1,1) ,KHFLIG (1,0, 1) , KHFLIG ( 1,1,1), 
KOFLIG (0,0,1),K0FLIG(0,1,1),K0FLIG(1,0,1),K0FLIG(1, 1,1) 



Same as above but for the diabatic blocking of proton single-particle parity-signature configu- 
rations. 
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2.8. Miscellaneous parameters 

This section lists keywords pertaining to definitions of parameters specifying various numerical 
conditions of the code execution. 



2.8.1. Keyword: MAX_MULTIP 

2, 4, 4 = NMUCON, NMUCOU, NMUPRI 

Maximum multipolarities A of multipole moments used in the code for the constraints, Eq. ()15|) . 
surface term of the Coulomb field, and printed on the output, respectively. Values not larger 
than A=9 are currently allowed. In case of the conserved parity, only even multipoles are used 
in the Coulomb field. 



2.8.2. Keyword: MAX_SURFAC 

0, = NSICON, NSIPRI 

Maximum multipolarities A of surface multipole moments, used in the code for the constraints, 
Eq. (fT3j) . and printed on the output, respectively. Values not larger than A=7 are currently 
allowed. 



2.8.3. Keyword: MAX_MAGNET 

0, = NMACON, NMAPRI 

Same as above but for the magnetic moments. In version (v2.08k) constraints on magnetic 
moments are not yet implemented, which requires NMACON=0. Values not larger than A=9 are 
currently allowed. 



2.8.4. Keyword: COULOMB 

80, 79, 0.25 = NUMCOU, NUMETA, FURMAX 

NUMCOU gives the number of points ]\['^°'^^, Eq. (I-lOO), used when summing up the Coulomb 
Green function. The dimensionless parameter d defining the size of the Coulomb parallelepiped, 
Eqs. (1-102) and (1-96), is given by NUMCOU* FURMAX. NUMETA gives the order of the Simpson 
integration of the solid harmonics on the faces of the parallelepiped, see Sec. 1-5. 



2.8.5. Keyword: SLOW-DOWN 

0.5, 0.5 = SLOWEV, SLOWOD 

The standard prescription to calculate the HF potential in the next iteration is to mix a given 
fraction e of the HF potentials from the previous iteration, with the fraction 1— e of potentials 
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calculated in the current iteration. SLOWEV and SLOWOD give the values of e separately for the 
time-even and time-odd potentials. This mixing slows down the convergence, and is necessary 
for the convergence of the iterative procedure. 



2.8.6. Keyword: SLOW-PAIR 

0.5 = SLOWPA 

Similarly as in the particle-hole channel, SLOWPA gives the value of the mixing fraction e used 
in the pairing channel, in analogy to SLOWEV and SLOWOD. 



2.8.7. Keyword: EPS_HERMIT 

l.OE-14 = EPSHER 

Numerical precision requested for determining values of Hermite polynomials. 



2.8.8. Keyword: OPTI_GAUSS 

1 = lOPTGS 

For I0PTGS=1 and IREAWS=0, expression = 2N^ + 2 is used to calculate the order of the 
Gauss-Hermite integration from the maximum number of HO quanta in the direction 
/i=a;, y, or z, and the input parameters NXHERM, NYHERM, and NZHERM are ignored. 



2.8.9. Keyword: GAUSHERMIT 

18, 18, 32 = NXHERM, NYHERM, NZHERM 

Orders L^, Ly, and Lz of the Gauss-Hermite integration in three Cartesian directions. Must 
be even. Ignored if IOPTGS=0 or IREAWS=1. 
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2.9. Parameters of the HO basis 

This section lists keywords pertaining to definitions of parameters of the HO basis. 
2.9.1. Keyword: BASIS_SIZE 

15, 301, 800.0 = NOSCIL, NLIMIT, ENECUT 

The HO basis is composed of states having not more than A^o=NOSCIL quanta in either of the 
Cartesian directions, and not more than M=NLIMIT states in total. The states are added to 
the basis according to the increasing energy of the deformed harmonic oscillator, 

^n.^nyn, = Tl^xiUx + \) + hujyiriy + |) + huj^{n^ + i). (12) 

In case of degenerate HO states (e.g., for an axially deformed HO) the complete multiplets 
are included, so the actual number of states LDBASE can be slightly larger than NLIMIT. If 
NLIMIT>0, the cutoff energy ENECUT is ignored; otherwise all states having HO energy smaller 
then ENECUT are included in the basis, and NLIMIT is ignored. 



2.9.2. Keyword: HOMEGAZERO 

1.2 = FCHOMO 

The code uses the standard value of the spherical HO frequency huo=4:lMeV / A^^^ multiplied 
by the scaling factor /=FCHOMO. 



2.9.3. Keyword: SURFAC_PAR 



86, 66, 1.23 



INNUMB, IZNUMB, ROPARM 



The code hfodd calculates parameters of the HO basis, and the zero-iteration Nilsson potential, 
by defining the standard nuclear surface S, 



S : R{9, 



c(a] 



1+ E E «A,nM(^,- 

A=0 fi.= ~\ 



(13) 



corresponding to the volume |7r*R0PARM^*(lNNUMB+IZNUMB), where c(a) is a function of a\fj, 
such that the volume enclosed by the surface S does not depend on a, i.e., in particular 
c(0) =i?o=ROPARM* ( I NNUMB+ IZNUMB) . 



2.9.4. Keyword: SURFAC_DEF 

2, 0, 0.61 = LAMBDA, MIU, ALPHAR 

The code defines frequencies of the deformed HO in three directions by using relation UxRx = 
ujyRy = uJzRz, where R^ are the lengths of principal axes of the nuclear surface E ([1^]) defined 
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by real deformation parameters a;v^=ALPHAR(LAMBDA,MIU). Tlie overall factor is defined by 
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2.10. Constraints on angular momentum 

This section lists keywords pertaining to definitions of constraints on angular momentum. The 
complete cranking term in the energy density functional is given by: 

/ c.,ox(Jo) V 



£ 



cran ^ ^ 



—^JOaiJoa) + Cja ({Joa) — Joa) — ^.Jla{Jla) 



2.10.1. Keyword: OMEGAY 

0.5 = OMEGAY 

Isoscalar angular frequency ujj^y in the y direction, see Eq. (fT^ . A non-zero value of OMEGAY 
requires IR0TAT=1 and cannot be used for the BCS pairing, i.e., for IPAIRI=1 and IPAHFB=0. 



2.10.2. Keyword: OMISOY 

0.00 = OMISOY 

Same as above for the isovector angular frequency ujjiy in the y direction. 



2.10.3. Keyword: OMEGA_XYZ 

0.00, 0.00, 0.00, = OMEGAX, OMEGAY, OMEGAZ, ITILAX 

For ITILAX=1, values of the three Cartesian components of the isoscalar angular frequency 
vector cjjo, see Eq. fll4j) . For ITILAX=0 these values are ignored. ITILAX=1 requires IR0TAT=1, 
ISIMPY=0, and IPAIRI=0. A non-zero value of OMEGAX, OMEGAY, or OMEGAZ requires broken 
symmetry ISIMTX=0, ISIMTY=0, or ISIMTZ=0, respectively. 



2.10.4. Keyword: OMEGA_RTP 

0.00, 0.00, 0.00, = OMERAD, OMETHE, OMEPHI, ITILAX 

Same as above but for values of the standard spherical components of the isoscalar angular 
frequency, i.e., wjor, c^joe, and lojq^. 



2.10.5. Keyword: OMISO_XYZ 

0.00, 0.00, 0.00, = OMISOX, OMISOY, OMISOZ, ITISAX 

For ITISAX=1, values of the three Cartesian components of the isovector angular frequency 
vector uJji, see Eq. (HH). For ITISAX=0 these values are ignored. ITISAX=1 requires IR0TAT=1, 
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ISIMPY=0, and IPAIRI=0. A non-zero value of OMISOX, OMISOY, or OMISOZ requires broken 
symmetry ISIMTX=0, ISIMTY=0, or ISIMTZ=0, respectively. 



2.10.6. Keyword: OMEGA_TURN 

= IMOVAX 

For IM0VAX=1 or —1, the isoscalar angular frequency vector ujjq is in each iteration set in the 
direction of the total angular momentum Jq, or opposite to this direction, respectively, while 
its length is kept equal to that fixed by the sum of squares of OMEGAX, OMEGAY, and OMEGAZ. For 
IMOVAX=0, vector ujjq is not changed during the iteration. IM0VAX=1 or —1 requires ITILAX=1, 
ITISAX=0, and IFLAGA=0. 



2.10.7. Keyword: SPINCONSTR 

0.0, 0.0, = STIFFI(2), ASKEDI(2), IFLAGI(2) 

For IFLAGI(2) = 1, the quadratic constraint on spin is used together with the linear constraint. 
Values of STIFFI(2) and ASKEDI(2) correspond respectively to Cjy and Joy in Eq. dUl). For 
IFLAGI(2)=0, there is no quadratic constraint on spin. This keyword is maintained only 
for compatibility reasons; in version (v2.08k) of the code hfodd it has been replaced by 
SPICON_XYZ. 



2.10.8. Keyword: SPICON_XYZ 

0.0, 0.0, = STIFFI(1),ASKEDI(1),IFLAGI(1) 

0.0, 0.0, = STIFFI (2) ,ASKEDI (2) ,IFLAGI (2) 

0.0, 0.0, = STIFFI (3) ,ASKEDI (3) ,IFLAGI (3) 

Same as above but for quadratic constraints in the x, y, and z directions corresponding to 
IFLAGI(1) = 1, IFLAGI(2)=1, and IFLAGI(1)=1, respectively. IFLAGI(1)=1, IFLAGI(2) = 1, 
or IFLAGI(3)=1 requires broken symmetry ISIMTX=0, ISIMTY=0, or ISIMTZ=0, respectively. 



2.10.9. Keyword: SPICON_OME 

0.0, 0.0, = STIFFA, ASKEDA, IFLAGA 

For IFLAGA=1, the quadratic constraint on the angle between the angular frequency and angular 
momentum vectors, see Eq. ()14|) . In version (v2.08k) the angle is constrained to zero. Value 
of STIFFA corresponds to Ca in Eq. (|T4|) . Value of ASKEDA must be set to 0; this variable is 
reserved for a future implementation of the constraint to a non-zero angle. For IFLAGA=0, there 
is no quadratic constraint on the angle. IFLAGA=1 requires ISIMPY=0, IR0TAT=1, ITILAX=1, 
ITISAX=0, IFLAGI(1)=0, IFLAGI(2)=0, IFLAGI(3)=0, and a non-zero value of sum of squares 
of OMEGAX, OMEGAY, and OMEGAZ. 
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2.11. Constraints on multipole moments 

This section lists keywords pertaining to definitions of constraints on multipole moments and 
surface multipole moments. The complete multipole constraint term in the energy density 
functional is given by: 

The multipole and surface multipole moments are defined as 

QA/i(r)=a,^rX(^>0). (16) 

and 

QUr) = a,,r'^'Y*^{e,<j,), (17) 

respectively, where Ya/i are the standard spherical harmonics in the convention of Ref. [Zj. For 
A<2, factors are listed in Table while for A>2 0^^4=1 . 

Table 5: Adopted definitions of the normalization factors and the corresponding multipole 
moments Qxfj,. 



X ji 


O-A/i 


Qa/1 







1 


1 


\/4vr/3 


z 


1 1 




X — iy 


2 


yi67r/5 


2^2 _ a;2 _ ^2 


2 1 


-^Stt/IS 


— zzy 


2 2 


y327r/5 


V3(x^ — — 2ixy) 



2.11.1. Keyword: MULTCONSTR 

2, 0, 0.01, 42.0, 1 = LAMBDA, MIU, STIFFQ, QASKED, IFLAGQ 

For IFLAGQ=1, the total multipole moment of the given multipolarity A and /i is constrained. 
Values of LAMBDA, MIU, STIFFQ, and QASKED correspond respectively to A, /i, Ca/j, and Qa/i in 
Eq. p5|) . For IFLAGQ=0, there is no constraint in the given multipolarity. 



2.11.2. Keyword: SURFCONSTR 

2, 0, 0.0, 0.0, = LAMBDA, MIU, STIFFS, SASKED, IFLAGS 
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For IFLAGS=1, the total surface multipole moment of the given multipolarity A and /i is con- 
strained. Values of LAMBDA, MIU, STIFFS, and SASKED correspond respectively to A, /i, Cf^, and 
Ql in Eq. For IFLAGS=0, there is no constraint in the given multipolarity. 
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2.12. Output-file parameters 

This section lists keywords pertaining to definitions of parameters specifying character of in- 
formation printed on the output file. 

2.12.1. Keyword: PRINT-ITER 

1, 0, 1 = IPRSTA, IPRMID, IPRSTO 

The code prints results for the first, middle, and/or last iteration if the corresponding param- 
eters IPRSTA, IPRMID, and IPRSTO equal 1. 



2.12.2. Keyword: PRINT-MOME 

1, 1, 1 = IPRI_N, IPRI_P, IPRI_T 

The code prints values of the neutron, proton, and total moments (multipole moments, surface 
multipole moments, or magnetic moments) only if the corresponding parameters IPRI_N=1, 
IPRI_P=1, and IPRI_T=1, respectively. 



2.12.3. Keyword: PRINT-INTR 

1 = INTRIP 

The code prints the values of moments (multipole moments, surface multipole moments, or 
magnetic moments) and angular momenta in the intrinsic frame only if INTRIP=1. 



2.12.4. Keyword: EALLMINMAX 

-12.0, 0.0 = EMINAL, EMAXAL 

The code prints tables of single-particle properties for states with values of the single-particle 
Routhians between EMINAL and EMAXAL. No table is printed unless EMINAL < EMAXAL. 



2.12.5. Keyword: EqUASI_MAX 

10.0 = EMAXQU 

The code prints tables of quasiparticle properties for the HFB states with values of the quasi- 
particle energies smaller than EMAXQU. No table is printed unless EMAXQU>0. 
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2.13. Files 

This section lists keywords that define files used by the code. 
2.13.1. Keyword: REVIEWFILE 

HFODD.REV = FILREV 

CHARACTER*68 file name of the review file. Must start at the 13-th column of the data line. 
The ASCII review file is written after calculating every data set (i.e., once per every EXECUTE 
item), provided IREVIE>0. The file is not rewound, so the results for several data sets can be 
accumulated in a single file. This is so provided the filename FILREV is not changed between 
the EXECUTE items. The file contains all relevant parameters and results of calculation in a 
form suitable for reading by another program. It is meant as an interface to programs which 
analyze and/or plot the results. The file contains sections defined by keywords (different than 
keywords used in the input data file described here). The detailed structure of the review file 
is not documented in the present write up, and can be inferred from inspecting the specimen 
produced by the sample run, and from the subroutine review. 



2.13.2. Keyword: REVIEW 

2 = IREVIE 

The review file is not written if IREVIE=0. For IREVIE=2, the table of single-particle properties 
is included in the review file in addition to other results, which are written for IREVIE=1. 



2.13.3. Keyword: REPLAYFILE 

HFODD.REP = FILREP 

CHARACTER*68 file name of the replay file. Must start at the 13-th column of the data line. The 
binary replay file with the name defined in FILREP must exist if IC0NTI=1, and will be read. If 
the filenames FILREP and FILREC are identical, the replay file will be subsequently overwritten 
as a new record file. These feature is implemented to facilitate chaining of jobs which follow 
one another. 



2.13.4. Keyword: RECORDFILE 

HFODD.REC = FILREC 

CHARACTER*68 file name of the record file. Must start at the 13-th column of the data line. If 
IWRIRE=1, binary record file is written after each HF iteration. It contains complete information 
which allows restarting the iteration in another run of the code. To restart, one has to specify 
IC0NTI=1 and provide the name of the file by defining FILREP. In case of the computer crash, 
or upon a successful completion of the given input data set, the record file contains the results 
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of the last performed iteration. The file is always rewound before it is written, so the results 
for consecutive iterations do not pile up. 



2.13.5. Keyword: RECORDSAVE 

1 = IWRIRE 

For IWRIRE=1, the record file, is saved on disc after each iteration is completed. For IWRIRE=0, 
it is saved only once, after all iterations are completed. This option is useful on systems where 
the speed of writing large amounts of data to disc could hamper the system performance. For 
IWRIRE=— 1, the record file is never saved, and the run cannot be later continued. 



2.13.6. Keyword: WOODSAFILE 

WOODS.WFN = FILWOO 

CHARACTER*68 file name of the Woods-Saxon file. Must start at the 13-th column of the 
data line. This file is read provided IREAWS=1. The binary Woods-Saxon file constitutes an 
interface between the Woods-Saxon code and the code hfodd. It contains the Woods-Saxon 
wave functions and numerous other parameters which define the current calculation. The 
parameters read from the Woods-Saxon file overwrite the values provided in the input data file. 
Since the current version of the Woods-Saxon code will be published separately, the feature of 
starting the iteration from the Woods-Saxon results is not documented in the present write up. 



2.13.7. Keyword: COULOMFILE 

HFODD. COU = FILCOU 

CHARACTER*68 file name of the Coulomb file containing auxiliary data for a faster calculation 
of the Coulomb field. Must start at the 13-th column of the data line. Parameters NUMCOU, 
NUMETA, and FURMAX which define the calculations of the Coulomb direct potential (see Sees. 
1-5 and II-3.5) are usually kept unchanged for the whole series of calculations in one region of 
nuclei. Therefore, many Coulomb auxiliary results can be calculated only once, and stored in 
the file FILCOU. Handling of this file is determined by the input parameters ICOULI and ICOULO. 



2.13.8. Keyword: COULOMSAVE 

0,0 = ICOULI, ICOULO 

Input parameters ICOULI and ICOULO determine actions pertaining to reading and/or writing 
of the Coulomb file FILCOU, according to the following table: 
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ICOULI ICOULO Action 

neitlier read nor write the Coulomb file 

1 read, but do not write the Coulomb file 

1 do not read, but write the Coulomb file 

1 1 use automated handling of the Coulomb file 

The default values of ICOULI=ICOULO=0 ensure that the code hfodd (v2.08k) behaves as 
that in version (vl.60r). However, unless it is required by special circumstances, a use of the 
automated handling of the Coulomb file (IC0ULI=IC0UL0=1) is recommended. Within the 
automated mode, the code checks whether the Coulomb file exists, and whether it contains 
data which match the current values of the input parameters NUMCOU, NUMETA, and FURMAX. If 
this is the case, the code reads the data from the Coulomb file. If this is not the case, the 
code calculates the Coulomb auxiliary results and stores them in the Coulomb file FILCOU. In 
the automated mode, the user is informed by appropriate messages printed on the output file 
about what type of the action has been taken in the given run of the code. 
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2.14. Starting the iteration 

This section lists keywords that pertain to starting and restarting the code. 



2.14.1. Keyword: RESTART 

= ICONTI 

For IC0NTI=1, results stored in the replay file (written as a record file in a previous run) are 
used to start the iteration. The code assumes (without checking) that the parameters of the 
HO basis are the same as those used in the previous run. The replay file name should be 
provided by defining FILREP. If the previous run was done with IREAWS=1, the current run 
must also use IREAWS=1, and the same Woods-Saxon file must be provided. This so because 
the Woods-Saxon file contains not only the information about the starting potential (which is 
ignored for IC0NTI=1) but also defines the HO basis. 



2.14.2. Keyword: CONT_PAIRI 

= IPCONT 

For IPC0NT=1, results stored in the replay file are used to define the HFB pairing properties in 
the first iteration; otherwise values read in matrices FERINI and DELINI are used. When the 
HFB pairing correlations are taken into account (i.e., for IPAHFB=1), and if a smooth restart 
and continuation of iterations from previously stored results is required, value of IPC0NT=1 
must be used. IPC0NT=1 is incompatible with either of IPAHFB=0 or ICONTI=0. 



2.14.3. Keyword: CONT_OMEGA 

= lOCONT 

For I0C0NT=1, the angular momentum vector Jq stored in the replay file is used to define in the 
first iteration the direction of the angular frequency vector uj jo, see Eq. (fT^ : otherwise values 
read in OMEGAX, OMEGAY, and OMEGAZ are used (see keyword OMEGA_XYZ). I0CGNT=1 requires 
IM0VAX=1 or -1. 



2.14.4. Keyword: READ_WOODS 

= IREAWS 

For IREAWS=1, the results stored by the previously performed Woods-Saxon calculation are 
used to start the iteration and to define the HO basis. For IC0NTI=1 and IREAWS=1 the 
Woods-Saxon file must also be provided, and is used only to define the HO basis. 
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2.14.5. Keyword: NILSSONPAR 

0, -1.175, -0.247, -1.175, -0.352, 11.17, 11.17, 6.28 = 
NILDAT, CNILSN, DNILSN, CNILSP, DNILSP, 
HBANIX, HBANIY, HBANIZ 

For IREAWS=0 and ICONTI=0, the code starts the calculation from the Nilsson potential. If 
NILDAT=1, the Nilsson parameters C and D [Ref. |Hj, Eq. (2.89)] are for neutrons given by 
CNILSN and DNILSN and for protons by CNILSP and DNILSP, while frequencies of the deformed 
HO are defined by HBANIX, HBANIY, and HBANIZ. If NILDAT=0, the Nilsson parameters C and 
D are defined by Eq. (2.91) and Table 2.3 of Ref. jH], while frequencies of the deformed HO are 
the same as those of the HO basis defined under keyword SURFAC_DEF. In the latter case the 
Nilsson parameters read from the input file are ignored. 
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3. OUTPUT FILE 

The code writes results on standard FORTRAN output file. The output file begins with the 
information pertaining to the general parameters of the calculation, then gives information 
about the starting point of the iteration, provides the convergence report, and finally contains 
the results calculated at the last iteration. Structure of the output file strongly depends on the 
type of calculation that are requested in the given run; especially on the assumed conserved 
and broken symmetries. As a rule, only the information relevant for the specific run is printed. 
Descriptions and comments printed on the output file are in most cases self-explanatory. In 
this section we list only those explanations that are not explicitly included in the output file. 

• Section CLASSICAL NUCLEAR SURFACE lists deformation parameters used to define the 
nuclear surface from which the basis parameters are derived, see keyword OMEGA_XYZ. It 
also gives the basis parameters such as the OSCILLATOR FREQUENCIES : HBAROX, HBAROY, 
and HBAROZ corresponding to huj^., huy, and huz- 

• Section PHYSICAL CONSTANTS gives values of he in MeVfm (H_BARC), of hc/e"^ (HBC0E2), 
of the neutron and proton masses in MeV/c^ (XMASSN and XMASSP), of the kinetic -energy 
coefficient /2m before (HBMASS) and after the center-of-mass correction (HBMRPA), of 
the elementary charge squared in MeVfm (ECHAR2), and of the coefficient preceding 
the integral in the Coulomb exchange energy (COULEX). For details see the comments in 
the SETBAS subroutine. 

• Section OSCILLATOR LENGTHS, CONSTANTS, and FREQUENCIES gives the values of 1/6^, 6^, 
and huj^, respectively, which characterize the HO basis in three Cartesian directions, ^=x, 

• Section BASIS CUT-OFF CONTROL PARAMETERS gives maximum numbers of the HO quanta 
in three directions NXMAXX, NYMAXX, and NZMAXX corresponding to A^^, Ny, and N^, as well 
as the orders of the Gauss- Hermite quadratures NXHERM, NYHERM, and NZHERM, correspond- 
ing to Lx, Ly, and L^, see keyword OPTI_GAUSS. It also gives the number M of the HO 
states included in the basis as requested in the input file (NLIMIT) and as used in the 
calculation (LDBASE). 

• Section SHAPE OF THE OSCILLATOR-BASIS DIAMOND gives the numbers of the HO quanta 
in a given direction for fixed numbers of the HO quanta in both remaining directions. 
The output is arranged in such a way that the shape of the grid of points n^nyUz is clearly 
visualized by projections in every of the three directions. 

• Section SKYRME FORCE DEFINITION lists the name and parameters of the Skyrme force, 
together with parameters KETA_J, KETA_W, KETACM, and KETA_M that define the way the 
given force should be used (see keyword SKYRME-STD). 

• Section COEFFICIENTS DEFINING THE SKYRME FUNCTIONAL gives values of the coupling 
constants in the Skyrme functional. These values take into account scaling factors, which 
are printed in the section SCALING FACTORS, unless they are all equal to 1. 
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• Section CALCULATIONS WITH THE TILTED-AXIS CRANKING gives values of components of 
the angular frequency vector and its length. For switches IM0VAX=1 or —1 and IOCONT=0, 
these values are used only in the first iteration, and later ignored, because the angular 
frequency vector is in each iteration readjusted to be aligned or anti-aligned with the 
angular momentum vector. 

• Sections PARITY/SIGNATURE, SIMPLEX, or PARITY CONFIGURATIONS give the vacuum and 
particle-hole configurations requested in the input file for the corresponding symmetry. 

• Section CONVERGENCE REPORT gives the list of performed iterations. For each iteration, 
one line is printed showing the energy, stability energy, average values of the intrinsic 
quadrupole moment Q and 7 deformation, total angular momentum, angular frequency, 
angle between the angular frequency and angular momentum vectors, and total pairing 
energy. 

• Section SINGLE-PARTICLE PROPERTIES lists the single-particle states calculated for the 
Nilsson, Woods-Saxon, or Hartree-Fock Routhian operators. For every state one line is 
printed which gives the value of the single-particle Routhian, the consecutive numbers 
in the parity-signature block or simplex block, the quantum numbers [A^, n2,A]f2 of the 
Nilsson asymptotic state which has the largest component in the given state, the average 
value of the parity operator (in %), the average values of projections of the intrinsic 
and total angular momenta (in Ti) , and their ratio called the (7- factor. For broken simplex 
symmetry, three projections of the total angular momentum and intrinsic spin are printed 
in the first and second line, respectively, for each single-particle state. 

• Section DENSITY INTEGRALS IN THE SKYRME FUNCTIONAL gives integrals of products of 
densities, which appear in the Skyrme functional. Terms in the functional are identified 
by the acronyms described under keywords EVE_SCA_TS and ODD_SCA_TS. 

• Section CONTRIBUTIONS TO ENERGY IN THE SKYRME FUNCTIONAL gives the values of var- 
ious terms which appear in the Skyrme functional. These contributions are the products 
of the coupling constants and of the density integrals described above. The sums of 
time-even and time-odd contributions are also printed. 

• Section EULER ANGLES OF THE PRINCIPAL-AXES FRAME gives the standard Euler angles 
a, /3, and 7 [H in degrees, which define the orientation of the intrinsic frame (principal- 
axes frame) of reference with respect to the original frame. 

• Section MULTIPOLE MOMENTS gives values of neutron, proton, or mass multipole moments 
with respect to the original frame of reference. In particular, the A=0 moment corre- 
sponds to the number of particles. Similarly, section MULTIPOLE MOMENTS [UNITS: (10 
FERMI) "LAMBDA] [INTRINSIC FRAME] gives analogous values with respect to the intrinsic 
frame of reference. Whenever the parity symmetry is broken, yet another analogous sec- 
tion gives information on values of multipole moments with respect to the center-of-mass 
reference frame. 

• Sections ROOT-MEAN-SQUARE AND GEOMETRIC SIZES gives the rms average values of the 
radius and of the x, and z coordinates. In order to better visualize the size of the 
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nucleus, the geometric sizes are also calculated by multiplying the rms radius by yS/S 
and the rms coordinates by y^- 

• Sections MAGNETIC MOMENTS [MAGNETON*FERMI" (LAMBDA-1)] and MAGNETIC MOMENTS 
[MAGNET0N*FERMI"(LAMBDA-1)] [INTRINSIC FRAME] give values of magnetic moments 
with respect to the original and intrinsic frame, respectively. 

• Section ANGULAR MOMENTA gives for ISIMPY=1 the average values of the total and intrin- 
sic neutron, proton, and total angular momentum (in h). It also gives the correspond- 
ing values and contributions to the first moment of inertia J^^^=I joj (in /i^/MeV). For 
ISIMPY=0, this section gives values of neutron, proton, and total projections of the total 
angular momentum and intrinsic spin on the three Cartesian axes in the original and 
intrinsic frame, respectively. 

• Sections NEUTRON CONFIGURATIONS and PROTON CONFIGURATIONS give a visual represen- 
tation of states occupied in the parity-signature blocks, simplex blocks, or parity blocks. 
The lines denoted by CONF: give the configurations requested in the input data set, while 
those denoted by VACC: give the configurations characterizing the given HF state. By 
comparing the two sets one can verify whether the requested configuration has been 
obtained, and eventually devise a new configuration to be calculated. The consecutive 
numbers printed in the horizontal direction correspond to the consecutive numbers in 
blocks printed in the section SINGLE-PARTICLE PROPERTIES. By comparing the two sec- 
tions one can effectively associate the Nilsson labels with the calculated configurations 
and also prepare the configuration input data. 

• Section ENERGIES gives a summary of the energies calculated for the HF state. The kinetic 
energy, single-particle energy, and pairing energy are printed for neutrons (NEU), protons 
(PRO), and all particles (TOT). Then, pairing rearrangement energy P-REARR, pairing 
gap PAIRGAP, and Fermi energy E-FERMI are printed. The Coulomb direct energy (DIR) 
and Coulomb exchange energy (EXC) are printed together with their sum (TOT). The 
multipole (MULT), surface multipole (SURF), and cranking (SPIN) constraint energies, 
are printed together with the corresponding corrections energies (CORR.). Then the 
rearrangement energy is printed followed by the values of the Routhian and the spin-orbit 
and Skyrme energies, the latter two split in the time-even (EVE) and time-odd (ODD) 
contributions. Finally, the total energies are printed as (SP) and (FUN), respectively, 
while their difference is printed as the stability (STAB). 
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FORTRAN source of code HFODD (v2.08k) is provided in the file hf odd . f and can be modified 
in several places, as described in this section. 

4.1. Dimensions of arrays 

The code HFODD uses the arrays' dimensions declared through the PARAMETER statements. This 
allows changing the dimensions and adapting the size of the reserved memory to the problem 
being solved. Whenever too small a dimension is defined, the code aborts with a message 
indicating the dimension that should be increased. Substantial amount of memory is required 
only for arrays depending on the following PARAMETER values: 

PARAMETER (NDMAIN=21) 

Corresponds to the maximum number of the HO quanta. Should be larger or equal to the input 
parameter NOSCIL defined under keyword BASIS_SIZE. 

PARAMETER (NDBASE=507) 

Corresponds to the maximum number of the HO basis states. Should be larger or equal to the 
input parameter NLIMIT defined under keyword BASIS_SIZE. It should also be larger or equal to 
the actual size of the HO basis LDBASE, which can be larger than NLIMIT in case of degenerate 
HO states. 

PARAMETER (NDXHRM=33 , NDYHRM=33 , NDZHRM=45) 

Correspond to the maximum numbers of the HO Gauss-Hermite integration nodes in three 
Cartesian directions. Should be larger or equal to the input parameters NXHERM, NYHERM, and 
NZHERM, respectively, that are defined under keyword GAUSHERMIT or determined in an optimal 
way, see keyword OPTI_GAUSS. 

PARAMETER (NDSTA=181) 

Corresponds to the maximum numbers of the single-particle states or quasiparticle states, which 
are kept after diagonalization. Should be larger or equal to the input parameters NUMBSP(0,0) 
+ NUMBSPd.O) and NUMBSP(0,1) + NUMBSP(1,1), see keyword PHASESPACE. 

On vector machines, parameters NDBASE, NDXHRM, NDYHRM, NDZHRM, NDSTAT and should be 
odd integers in order to minimize the risk of bank memory conflicts. By the same token, 
parameter NDMAIN should be even, because it defines the matrix dimensions beginning with 0. 

4.2. Vectorization properties 

As discussed in Sec. 1-4.3, the code HFODD has to operate by using seven-fold nested short 
loops, and this part does not perform well in a vector processor. However, it turns out that 
the loops can be artificially made longer in such a way that the final CPU time in a vector 
processor actually becomes much shorter. All the places where this trick has been applied can 
be identified in the source file by finding the lines beginning with CVECTOR and CSCALAR, for 
example: 

CVECTOR 

DO KZ=0 , LAZOXY (NX , NY) +LAZOXY (MX , MY) 
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CSCALAR DO KZ=NZ+MZ , , -2 

RESULT=RESULT+COEFOO (KZ , NZ ,MZ , 3) *T_AUXI (KZ) 
END DO 

The line beginning with CSCALAR should be made active on a scalar or superscalar machine. 
The line immediately below the line which begins with CVECTOR should be active on a vector 
machine. The results of calculation do not depend on which version of loops is activated. 

4.3. Library subroutines 

The code hfodd requires an external subroutine which diagonalizes complex hermitian matri- 
ces. Version (vl.60r), see II, has been prepared with an interface to the NAGLIB subroutine 
f02axe, and version (vl.75r), see III, with an interface to the LAPACK subroutine zhpev. 
In the present version (v2.08k), both these interfaces remain supported, and can be activated 
as described in II and III, respectively. However, the recommended interface is now to the 
LAPACK subroutine zhpevx, as described in this section. 

In version (v2.08k) we have implemented interface to the LAPACK subroutine zhpevx, 
which can be downloaded (with dependencies) from 

http://netlib2.cs.utk.edu/cgi-bin/netlibfiles.pl?filen"ame=/lapack/complexl6/zhpevx.f| 
This subroutine finds not all, but only the lowest eigenvectors, and hence performs calcula- 
tions in less CPU time; the gain is particularly significant for a large HO basis. Numbers of 
eigenvectors to be found are defined by the size of the phasespace, see keyword PHASESPACE. 

Subroutine zhpevx makes use of the BLAS (Basic Linear Algebra Subprograms) library, 
which can be included in the above downloaded file (unoptimized version). A version optimized 
for a given computer can be used as well, if it is available in the system. 

Subroutine zhpevx and its dependencies are in the REAL*8/C0MPLEX*16 version, and should 
be compiled without promoting real numbers to the double precision. On the other hand, the 
code hfodd itself does require compilation with an option promoting to double precision. 
Therefore, the code and the zhpevx package should be compiled separately, and then should 
be linked together. 

In order to activate the interface to the LAPACK zhpevx subroutine, the following modi- 
fications of the code hfodd (v2.08k) have to be made: 

1. Change everywhere the value of parameter I_CRAY=1 into I_CRAY=0. 

2. Change everywhere the value of parameter IZHPEV=0 into IZHPEV=2. 

3. If your compiler does not support undefined externals, or subroutines called with different 
parameters, remove calls to subroutines cgemm, f02axe, and zhpev. 

A set of c-shell and ex-editor scripts is provided within the hfodd distribution file, which 
allows an easy installation, compilation, and execution of the code on a Linux computer running 
the INTEL© FORTRAN COMPILER. 



4.5. Printing the execution times 
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4.4. FORTRAN-90 version 

The code hfodd version (v2.08k) is written in FORTRAN- 77. However, on several platforms 
several tested FORTRAN-77 compilers did not perform well with respect to the memory man- 
agement. In particular, they often distributed large matrices on computer's stack and heap in 
a wasteful way, and then the available memory could be insufficient for the code to execute, 
especially for a large HO basis. Therefore, crucial parts of the code have been promoted to 
FORTRAN-90, with a use of an explicit memory allocation and deallocation. In the FORTRAN 
source code provided in the file hfodd. f, all these FORTRAN-90 features are commented out 
and inactive. However, very simple modifications of the source code can easily be performed 
to transform code HFODD to FORTRAN-90. To this end, the user should run an automated 
editor script which performs the following two operations: 

1. Replace all strings "C F90" by spaces, and 

2. Remove two lines of code (or three lines of the file) after lines marked by the string: 
"C IN F90 VERSION PLEASE REMOVE THE TWO FOLLOWING LINES". 

The FORTRAN-90 version should be compiled with the "fixed form" non-standard option of 
the FORTRAN-90 compiler. 

A set of c-shell and ex-editor scripts is provided within the hfodd distribution file, and 
allows performing the above modifications of the FORTRAN source code. 

4.5. Printing the execution times 

After successful completion of execution, the code hfodd prints the table of CPU times spent 
in principal subroutines, as well as the total execution time. Since implementation of the 
time-checking functions widely varies between compilers and platforms, in the file hfodd. f all 
references to these functions are commented out and inactive. However, very simple modifica- 
tions of the source code can easily be performed to activate one of the following options: 

1. Replace all strings "CCPUT" by spaces in order to activate a call to the FORTRAN-90 
subroutine cpu_time, which returns CPU execution times. 

2. Replace all strings "CETIM" by spaces in order to activate a call to function etime, 
which returns CPU execution times. 

3. Replace all strings "CUNIX" by spaces in order to activate a call to the UNIX function 
TIME, which returns wall clock execution times. 

4. Replace all strings "CSAL" by spaces in order to activate a call to the SALFORD© 
FORTRAN subroutine dclock(3), which returns wall clock execution times. 

5. Replace all strings "CF32" by spaces in order to activate a call to the MICROSOFT© 
F32 subroutine gettim, which returns wall clock execution times. 

Other time-checking functions, available to the user, can be analogously called in the subroutine 

CPUTIM. 
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5. TEST RUNS 

Several examples of input and output files are provided with the released versions, (vl.60r), 
(vl.75r), and (v2.08i), of the program HFODD, or added in this user's guide, as listed in Table IHl 
The input data files dyl52-a.dat through dyl52-e.dat present the valid input items of each 
version and the default or recommended values of the input parameters. The input data files 
tbl51-a.dat through tbl51-d.dat are distributed with version (vl.75r), and provide examples 
of calculations with diabatic blocking. The input data files lal32-a.dat and lal32-b.dat 
are specific for version (v2.08i), and illustrate various symmetry options of the code. The files 
pbl28-a.dat and snl20-a.dat are added in this user's guide and give examples of the spherical 
HF and HFB calculations, respectively. The test runs are described in detail in Section 15.11 
and the sample input data files are reprinted in Section EIH 



Table 6: Input data files supplied with versions (vl.60r), (vl.75r), and (v2.08i) of the program 
HFODD or added in this user's guide. Each output file listed has been obtained by running 
version (2.08k) of the code for corresponding input data file. 



Distribution 


Input 


Contents 


Output 


(vl.eOr) 


dyl52- 


-a 


dat 


Only keyword EXECUTE 

Example of superdeformed HF calculations 


dyl52-a 


out 




dyl52-b.dat 


All input items of (vl.60r). 


dyl52-b 


out 










default values of input items 






(vl.75r) 


dyl52- 


-c 


dat 


All input items of (vl.75r), 
default values of input items 


dyl52-c 


out 




dyl52- 


-d 


dat 


Input items added in (vl.75r), 
recommended values of input items 


dyl52-d 


out 




tbl51- 


-a 


dat 


Example of adiabatic calculations 


tbl51-a 


out 




tbl51- 


-b 


dat 


Example of adiabatic calculations 


tbl51-b 


out 




tbl51- 


-c 


dat 


Example of diabatic blocking 


tbl51-c 


out 




tbl51- 


-d 


dat 


Example of diabatic blocking 


tbl51-d 


out 


(v2.08i) 


dyl52- 


-e 


dat 


All input items of (v2.08i), 
default values of input items 


dyl52-e 


out 




lal32- 


-a 


dat 


Example of tilted-axis-cranking calculations 


lal32-a 


out 




lal32- 


-b 


dat 


Example of tilted-axis-cranking calculations 


lal32-b 


out 


User's 


pbl28- 


-a 


dat 


Example of spherical HF calculations 


pb208-a 


out 


guide 


snl20- 


-a 


dat 


Example of spherical HFB calculations 


snl20-a 


out 



5.1. Descriptions of the test runs 
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5.1. Descriptions of the test runs 
5.1.1. Files dyl52-a.dat through dyl52-e.dat 

As explained in Section all the input data parameters of the code have their default values, 
which are assumed if the relevant input items are absent in the input data file. Each parameter 
introduced in a given version is also valid in all later versions, and its default value remains 
unchanged. Later versions introduce new parameters, but their default values are such that the 
new features are not active. Thus, the results of running a later version are identical to those 
of running an earlier version, provided that the values of those new parameters are left at the 
default. In particular, running each version with its default data set corresponds to identical 
conditions and gives the same solution. All the input items valid for a given version are also 
recognized and read by the subsequent ones. In general, running a later version with an input 
data file prepared for an earlier version leads to the same solution. However, the structure of 
the output file may be modified from version to version; in particular because later versions 
may calculate more observables. 

Running the code with the default values results in performing 50 iterations for the super- 
deformed state in ^^^Dy at the angular frequency of hu = 0.5 MeV. The pairing correlations 
are not taken into account. The HO basis is defined by an axially symmetric (with respect 
to the z axis) nuclear surface with physical deformation typical for superdeformed bands in 
the A ~ 150 region. The parity, y-signature, x-simplex^ (and their products) are imposed as 
conserved symmetries, while the time-reversal symmetry is broken in order to describe rota- 
tion (about the y axis). The run starts from the Nilsson potential, and an appropriate weak 
constraint on the quadrupole moment (Q20) is used to ensure that the iterations fall into the 
superdeformed minimum. 

The input data files dyl52-a.dat through dyl52-e.dat define either the default or recom- 
mended values of the input parameters for the released versions of the code. They contain either 
all or some of the valid input items, as described in Table El File dyl52-a.dat contains only 
keyword EXECUTE, and when used as input for any version of the code it results in performing 
calculations with the default data set. Files dyl52-b.dat, dyl52-c.dat, and dyl52e-dat con- 
tain all the valid input items of versions (vl.60r), (vl.75r), and (v2.08i), respectively, and the 
input data are identical to the default values. File dyl52-d.dat contains only the input items 
introduced in version (vl.75r), and the corresponding new parameters are set to recommended 
values, which differ from the default ones in the following way. The default values are such that 
the code works identically to version (vl.60r). The parameters EPSITE, NULAST, EPSPNG, and 
NUCHAO, control the termination of the iterations. By default, they are set so that always the 
number of iterations specified by NOITER is performed. The recommended values define an op- 
timal level, at which convergence, "chaotic" divergence, or "ping-pong" divergence are detected 
to stop the iterations earlier. The differences between the default and recommended values of 
the variables governing the diabatic blocking are irrelevant because the diabatic blocking is 
disabled by default. Also by default, the Coulomb file is neither read from nor written to disc. 
With the recommended values of the parameters ICOULI and ICOULO, its reading and writing 
is handled automatically, which saves the CPU time. 

Files dy 152-a . dat through dyl52-e . dat serve only as suitable patterns to modify the input 
parameters. However, it should be a good practice to include in the input data file only those 
input items which modify the input parameters with respect to the default values. In this 
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J I I I I I I II I I I I I I L 



0.5 0.6 0.7 0.8 0.5 0.6 0.7 0.8 
Rotational frequency (MeV) 

Figure 1: Negative-parity neutron single-particle Routhians in ^^^Tb calculated for the 
[761]3/2(r=+i) (a) and [514]9/2(r=+i) (b) neutron diabatic configuration. Solid and dashed 
curves denote the r=+i and r=-i signatures, respectively. The arrows denote the angular 
frequencies where the converged solutions near the crossing points could not be found. 

way, the danger of accidentally modifying some intricate numerical parameters of the code is 
minimized. 

5.1.2. Files tbl51-a.dat through tbl51-d.dat 

Files tbl51-a.dat through tbl51-d.dat provide data for calculations in which single-particle 
configurations corresponding to two crossing bands in ^^^Tb jHI are followed either adiabatically 
or diabatically. The pairing correlations are not taken into account. The parity, y-signature, 
x-simplex^ (and their products) are imposed as conserved symmetries, while the time-reversal 
symmetry is broken in order to describe rotation (about the y axis). The angular frequency 
range from hu = 0.50 to 0.80 MeV is scanned with the step of 0.05 MeV. Calculations at 
hu = 0.50 are started from the Nilsson potential, and iterations for each subsequent frequency 
are restarted from the solution obtained at the preceding frequency. The two considered bands 
are based on excited physical configurations, in which the last neutron occupies either the 
[514]9/2(r=+i) or [761]3/2(r=+i) Nilsson level. At hoo = 0.50 MeV the level [514]9/2(r=+i) 
is lower in energy, and has number 23 --|- (23rd state of 7r=-l and r=+i), while the level 
[761]3/2(r=-|-i) has number 24--I-. These two levels cross in the concerned frequency interval 
so that at hu = 0.80 MeV the situation is the opposite. This is illustrated in Figure ^ which 
is explained in detail below. 



5.1. Descriptions of the test runs 
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In the files tbl51-a.dat and tbl51-b.dat it is required that always the level 23-+ or 24- 
+ be occupied, respectively, which amounts to following the configuration adiabatically. If, say, 
the level 23-+ is occupied, the solutions at low and high frequencies correspond to physical 
configurations [514]9/2(r=+z) and [761]3/2(r=+z), respectively, while in the vicinity of the 
crossing no solution can be found, because the so-called "ping-pong" divergence occurs. This 
means that the physical states [514]9/2(r=+z) and [761]3/2(r=+«) change their energy order 
every second iteration, so that each of them once becomes the level 23-+. Thus, putting the 
neutron always on the level 23-+ results in occupying alternately the very different physical 
states [514]9/2(r=+i) and [761]3/2(r=+2), which cannot lead to convergence. 

The code detects the "ping-pong" divergence, stops the iterations, and prints a table that 
helps identify the crossing orbitals. Such a table is reproduced in Table[7|for the 23- + adiabatic 
configuration at hu = 0.65 MeV. For both protons and neutrons, in each parity-signature block, 
and separately for the particle and hole states, the code finds the level for which the difference 
in the angular momentum or intrinsic spin alignment (on the y axis) between the last and last 
but one iteration is the largest. The indices of those levels in the parity-signature blocks, and 
the absolute values of the differences are printed in the upper part of the table, independently 
for the angular momentum and intrinsic spin. If in some symmetry block, the indices of the 
particle and hole states found in this way differ by one, then it is probable that these two 
levels cross and cause the "ping-pong" divergence. When such pairs exist, they are listed in the 
lower part of the table. In the given example, the pair of states 23-+ and 24-+ is properly 
recognized on the basis of both the angular momentum and intrinsic spin. The crossing orbitals, 
[514]9/2(r=+i) and [761]3/2(r=+z), have markedly different alignments of about —0.46^ and 
+2.23, respectively, which corresponds to the values 2.820 and 2.807/i given in the table. 

If one wishes to follow, say, the configuration [761]3/2(r=+i) diabatically, one should request 
that, of the levels 23-+ and 24-+, the one with physical properties of the state [761]3/2(r=+z) 
be occupied. One can chose, for instance, the state with higher angular momentum alignment, 
independent of whether it is lower or higher in energy. In order to diabatically follow the 
configuration [514]9/2(r=+z), one can chose the state with lower alignment. The corresponding 
runs are defined in the files tbl51-c.dat and tb-151-d.dat. Indeed, in each case solutions 
corresponding to the same physical configuration is found at each angular frequency. Moreover, 
the "ping-pong" divergence does not occur, even very close to the crossing frequency. Figure [T] 
shows the negative-parity neutron single-particle Routhians obtained from the calculations of 
files tbl51-c.dat and tb-151-d.dat. 

5.1.3. Files lal32-a.dat and lal32-b.dat 

File lal32-a.dat defines two consecutive runs, which lead to a planar solution fO] in the 
triaxially deformed nucleus ^'^^La for hu = 0.25 MeV. The pairing correlations are not taken 
into account. Since the calculations involve the tilted-axis cranking, a spherically invariant 
HO basis (with the three HO frequencies equal and comprising only entire HO shells) is used 
in order not to favor any orientation of the solution in space. The first run is performed at 
hu = with the parity, ^/-signature, x-simplex^ (and their products) conserved, and with the 
time-reversal broken in order to describe odd numbers of neutrons and protons. Due to the 
presence of the odd nucleons, even for huj=0 one obtains a non-zero angular momentum, which 
is oriented along the y axis. The run starts from the Nilsson potential, and in 94 iterations finds 
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Table 7: Example of the output printed when the "ping-pong" divergence is found, see text. 

* * 

* TABLE BELOW GIVES THE MAXIMUM ABSOLUTE VALUES OF DIFFERENCES OF ALIGNMENT * 

* BETWEEN THE LAST TWO ITERATIONS. IT MAY SERVE AS A GUIDE TO SELECT STATES * 

* WHICH EXCHANGE WAVE FUNCTIONS IN EVERY SECOND ITERATION. SUCH STATES MAY * 

* BE CANDIDATES FOR A CALCULATION WITH DIABATICALLY BLOCKED CONFIGURATIONS * 

* * 

* * * 

* BLOCKS * ANGULAR MOMENTUM ALIGNMENT INTRINSIC SPIN ALIGNMENT * 

* * PARTICLES HOLES PARTICLES HOLES * 

* ISO. PAR. SIG. * INDEX VALUE INDEX VALUE INDEX VALUE INDEX VALUE * 

* * * 



* 








* 


























* 


* 


NEUT 


+ 


+ 


* 


38 


0. 


.221 


13 


0. 


.018 


24 


0. 


.023 


13 


0. 


.016 


* 


* 


NEUT 


+ 




* 


28 


0. 


.163 


7 


0. 


.041 


33 


0. 


.014 


7 


0. 


.032 


* 


* 


NEUT 






* 


39 


0. 


.095 


13 


0. 


.198 


44 


0, 


.015 


14 


0. 


.150 


* 


* 


NEUT 




+ 


* 


24 


2, 


.820 


23 


2. 


.807 


24 


0. 


.164 


23 


0. 


.167 


* 


* 


PROT 


+ 


+ 


* 


24 


0. 


.255 


13 


0. 


.149 


25 


0. 


.129 


14 


0. 


.122 


* 


* 


PROT 


+ 




* 


27 


1. 


.733 


13 


0. 


.062 


21 


0. 


.319 


8 


0. 


.020 


* 


* 


PROT 






* 


32 


0. 


.119 


10 


0. 


.530 


32 


0. 


.029 


10 


0. 


.065 


* 


* 


PROT 




+ 


* 


24 


2. 


.649 


15 


0. 


.082 


24 


0. 


.153 


10 


0. 


.030 


* 


* 








* 


























* 



* * 

* FOLLOWS THE LIST OF CANDIDATE CONFIGURATIONS FOR THE DIABATIC BLOCKING * 

* * 

* NEUT - + * 24 2.820 23 2.807 * 

* NEUT - + * 24 0.164 23 0.167 * 

a triaxial solution, whose shape is forced by using appropriate constraints for the quadrupole 
moments (Q20) and {Q22)- The constraints are set so that this solution nearly corresponds to 
the unconstrained minimum, and are used only to prevent the iterations from falling into other 
local minima. These constraints are chosen so that the short, medium, and long axes of the 
mass distribution coincide with the y, x, and z axes, respectively. For the concerned solution, 
the single-particle configuration defined in the input data file in terms of numbers of particles 
in the parity-signature blocks physically corresponds to the odd proton particle occupying the 
lowest substate of the /iii/2 orbital and the odd neutron hole left on the highest /iii/2 substate. 

The second run starts from the solution obtained in the first run, and performs calculations 
with all symmetries broken except parity, with multipole constraints released, and for the 
angular frequency of hu = 0.25 MeV. The single-particle configuration specified in terms of 
parity blocks is obtained from the configuration for the first run simply by adding the numbers 
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of particles in tlie two signature sub-blocks of each parity block. All the three components of 
the initial angular frequency vector are non-zero, hux = hujy = huiz = (0.25/\/3) MeV, which 
breaks all symmetries except parity in the first iteration, thus making the calculations otherwise 
symmetry-unrestricted. If a symmetry was not broken in the first iteration, it would artificially 
remain conserved as governed by the theorem of self-consistent symmetries jH] . 

The Kerman-Onishi condition jTT| states that the angular momentum and angular frequency 
vectors must be parallel in self-consistent solutions. Standard iteration, with IMOVAX=0, re- 
quires several thousand iterations for the nucleus to turn appropriately in space to achieve 
the alignment of these vectors below the angle of 0.01°. On the other hand, iteration with 
IM0VAX=1, as shown in the example of file lal32-a.dat, achieves the alignment of 0.000001° 
within 92 iterations. In the solution resulting from the second run, the angular momentum 
vector has non-zero components only on the y and z axes of the intrinsic frame, which are the 
short and long axes of the triaxial mass distribution. The calculations of the file lal32-a.dat 
are performed in two steps to get a better control over what is going on and to avoid divergence 
or falling into undesired minima. 

The input data provided in the file lal32-b . dat are similar to those of the file lal32-a. dat, 
except that the second run is performed with the x-simplex"^ conserved. In this way, the 
angular momentum and angular frequency vectors are confined within the plane y-z, which 
is also a principal plane of the mass distribution. After the first run, this plane contains 
the short and long axes. The components of the initial angular frequency vector are set to 
hujy = hu, = (0.25/v/2)MeV and hux=0. Since the planar solution in question corresponds 
to non-equal projections of the angular momentum on the short and long axes, in order to 
satisfy the Kerman-Onishi condition, the nucleus still has to turn in space about the x axis. 
For IMOVAX=0, this process is also very slow, while for IM0VAX=1, as in file lal32-b.dat, 
it takes again only 92 iterations. The x-simplex-^-imposing solution obtained by using the 
file lal32-b.dat differs from the symmetry-unrestricted result of the file lal32-a.dat only by 
orientation in space; apart from that, the values of all observables are identical. This shows that 
in the considered self-consistent solution the x-simplex-^ symmetry is spontaneously conserved. 

5.1.4. File pb208-a.dat 

The input data file pb208-a.dat defines two consecutive runs to perform precise calculations 
for the spherical ground state of the doubly-magic nucleus ^osp]-, large basis of 26 HO 
shells. The pairing correlations are not taken into account. Conservation of all the symmetries 
considered in the code (the group -Dj^) is imposed. For the spherical shape, the specified 
parity-signature configuration corresponds to occupation of the entire magic shells, up to the 
proton and neutron magic numbers 82 and 126, respectively. A spherical basis is used, i.e., 
with all the three characteristic HO frequencies equal {huJx=fiujy=huJz=fiu!o) and comprising 
only the entire HO shells. The so-called physical value of hcuQ = / x 41 MeV/A^^^ with / = 1.2 
is taken, following the default settings of the code. To speed up the calculations, the program 
HFODD uses a diagonalization algorithm that can reduce the computational effort by finding 
less eigenstates than the matrix size (the lowest ones). In the present example, in each charge- 
simplex block only one extra state is calculated as compared to the number of occupied states in 
that block, as governed by the input parameter NUMBSP. This is far less than the code's default, 
and such a way of proceeding is possible only because no particle-hole excitations into higher 
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states are required. Also for the sake of saving time, a solution in a small basis of iVo=12 HO 
shells is obtained first, starting from the Nilsson potential, and then serves as a good starting 
point for the time-consuming calculations with Nq=26. The number of states in the spherical 
basis comprising the HO shells from to No equals M={Nq + l)(iVo + 2)(A^o + 3)/6, and hence 
the numbers of the HO basis states (for each projection of the intrinsic spin) equal M=455 and 
M=3654, respectively for iVo=12 and Nq=2Q. Since the code hfodd can restart iterations only 
if the number of Gauss-Hermite integration points is preserved, both runs are performed with 
54 points in each of the Cartesian directions, which is excessive for Nq=12, but optimal for 
Nq=2Q. For Nq=2Q, the 54 points ensure exact integration of all matrix elements apart from 
those of the Coulomb interaction and density-dependent term. 

The energies of ^osp^ obtained from calculations with 7Vo=12 and A'^o=26 equal —1630.638 
and —1635.531 MeV, respectively. The latter value is very close to the result of —1635.666 MeV 
obtained by using a very precise one-dimensional spherical code working in spatial coordinates 
|12j . Although the Nq=\2 HO basis is evidently too small for a precise determination of the 
energy, other observables may not be so sensitive to the HO basis size. For instance, the mass 
root-mean-square radii obtained for A''o=12 and A^o=26 equal 5.5540 and 5.5554 fm, respectively, 
the result of coordinate-space calculations being 5.5550 fm. 

5.1.5. File snl20-a.dat 

The input data file snl20-a.dat provides an example of the HFB calculations for the spherical 
ground state in ^^°Sn. All symmetries considered in the code (the group -D^^) are imposed, 
except the particle-number symmetry, which is broken in order to describe the pairing corre- 
lations. A spherical basis is used. The pairing correlations are included within the full HFB 
method, and the zero-range density- dependent interaction given by the form-factor (lllj) is used 
in the particle-particle channel. The particle and pairing densities are obtained by summing 
up contributions from the quasiparticle states with the equivalent-spectrum energies not ex- 
ceeding emax=60MeV, which is the code's default pairing window. In the form-factor pi|) . 
the power a=l and the saturation density po=0.32fm~^ are taken. The interaction strength, 
Vo=— 285.88 MeV fm^, was adjusted for the given cutoff energy, emax, to reproduce the value of 
Ajv=1.245MeV for the average neutron pairing gap in ^^°Sn. The same values of a, po and Vq 
are taken for both neutrons and protons. The Fermi energy of —8 MeV and the pairing gap of 
1 MeV are used in the first iteration to construct the HFB Hamiltonian, while the initial mean 
field is given by the Nilsson potential. 

Since ^^°Sn is a proton-magic nucleus, the proton pairing decreases during the iteration and 
vanishes in the HFB solution. As expected, one obtains a significant thickness of the neutron 
skin (difference between the neutron and proton root-mean-square radii), which in the present 
case equals 0.14 fm. 
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5.2. Listings of the test input data files 

5.2.1. File dyl52-a.dat 

EXECUTE 
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5.2.2. File dyl52-b.dat 



General data 



NUCLIDE 
ITERATIONS 

SKYRME-SET 
EVE_SCA_TS 
ODD_SCA_TS 
EVE_SCA_PM 
ODD_SCA_PM 
G_SCALING 

ROTATION 
SIGNATUREY 
TSIMPLEXES 
PAIRING 

VACSIM_NEU 
VACSIM_PRO 
PHSIMP_NEU 
PHSIMP_PRO 
VACSIG_NEU 
VACSIG.PRO 
PHSIGN_NEU 
PHSIGN_PRO 

MAX_MULTIP 



86 66 
50 



SKM* 

RHO 
1. 1. 

SPI 
1. 1. 

RHO 
1. 1. 

SPI 
1. 1. 

1.0 1.0 



RHOD 
1. 1. 

SPID 
1. 1. 

RHOD 
1. 1. 

SPID 
1. 1. 



1 
1 




Interaction 



LPR 
1. 1. 

LPS 
1. 1. 

LPR 
1. 1. 

LPS 
1. 1. 



TAU 
1. 1, 

CUR 
1. 1, 

TAU 
1. 1, 

CUR 
1. 1, 



SCU 
1. 1. 

KIS 
1. 1. 

SCU 
1. 1. 

KIS 
1. 1. 



Symmetries 



Configurations 



SIMP SIMM 

43 43 
SIMP SIMM 
33 33 

PART HOLE 
00 00 00 00 

PART HOLE 
00 00 00 00 

PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
16 16 17 17 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 00 00 



00 



00 



HOLES 
00 00 

HOLES 
00 00 



DIV 
1. 1, 

ROT 
1. 1, 

DIV 
1. 1, 

ROT 
1. 1. 



00 
00 



Numerical data 
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COULOMB 

SLOW-DOWN 

EPS_HERMIT 

OPTI_GAUSS 

GAUSHERMIT 

BASIS_SIZE 
HOMEGAZERO 
SURFAC_PAR 
SURFAC_DEF 
SURFAC_DEF 

OMEGAY 

MULTCONSTR 

SPINCDNSTR 

PRINT- ITER 
EALLMINMAX 

REVIEWFILE 

REVIEW 

RECDRDFILE 

REPLAYFILE 

WDDDSAFILE 

RESTART 



2 4 4 

80 79 0.25 
0.5 0.5 
l.OOE-14 
1 

18 18 32 



15 301 
1.2 

86 66 

2 

4 



0.50 
2 



— Parameters of the HO basis 



800. 

1.23 
0.61 
0.10 



Constraints 



0.01 42. 1 

0.00 0. 

Output-file parameters 



-12. 0. 

HFODD . REV 
2 

HFODD . REC 
HFODD . REP 
WOODS . WFN 



Files 



Starting the iteration 
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READ_WDDDS 



NILSSONPAR 

-1.175 -0.247 -1.175 -0.352 11.170 11.170 6.280 
Calculate 

EXECUTE 

Terminate 

ALL_DONE 



5.2. Listings of the test input data files 
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5.2.3. File dyl52-c.dat 



This file (dyl52-c.dat) contains the input data for the code HFODD. | 
Values of all parameters are identical to the default values assumed by | 
the code when no input data are provided. I 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

PHASESPACE 



86 66 

50 

0.0 

1 

0.0 3 




General data 



Interaction 



SKYRME-SET 

SKM* 

EVE_SCA_TS RHD 

1. 1. 
ODD_SCA_TS SPI 

1. 1. 
EVE_SCA_PM RHO 

1. 1. 
ODD_SCA_PM SPI 

1. 1. 

G.SCALING 

1.0 1.0 



ROTATION 
SIGNATUREY 
TSIMPLEXES 
PAIRING 

VACSIM_NEU 
VACSIM PRO 



1 
1 
1 




RHOD 
1. 1. 

SPID 
1. 1. 

RHOD 
1. 1. 

SPID 
1. 1. 



LPR 
1. 1. 

LPS 
1. 1. 

LPR 
1. 1. 

LPS 
1. 1. 



TAU SCU 

1. 1. 1. 1. 

CUR KIS 

1. 1. 1. 1. 

TAU SCU 



1. 1. 

CUR 
1. 1. 



Symmetries 



1. 1. 
KIS 
1. 1. 



Configurations 



SIMP SIMM 

43 43 
SIMP SIMM 



DIV 
1. 1, 

ROT 
1. 1, 

DIV 
1. 1, 

ROT 
1. 1, 
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PHSIMP.NEU 
PHSIMP_PRO 
VACSIG_NEU 
VACSIG_PRO 
PHSIGN.NEU 
PHSIGN_PRO 
DIASIM_NEU 
DIASIM_PRO 
DIASIG_NEU 
DIASIG.PRO 

MAX_MULTIP 

COULOMB 

SLOW-DOWN 

EPS_HERMIT 

OPTI .GAUSS 

GAUSHERMIT 

BASIS.SIZE 
HOMEGAZERO 
SURFAC_PAR 
SURFAC_DEF 
SURFAC_DEF 

OMEGAY 



33 33 
PARTICLES 

00 00 
PARTICLES 
00 00 
PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
16 16 17 17 
PARTICLES 

00 00 
PARTICLES 

00 00 
PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 
02 02 



00 



00 



02 



02 



00 



00 



02 



02 



80 79 0.25 

0.5 0.5 
l.OOE-14 
1 

18 18 32 

15 301 800. 
1.2 

86 66 1.23 

2 0.61 

4 0.10 



HOLES 

00 00 

HOLES 
00 00 



00 



00 



01 



01 



HOLES 
00 00 
HOLES 

00 00 
HOLES 

01 01 
HOLES 

01 01 
HOLES 

01 01 
HOLES 

01 01 



Numerical data 



00 
00 



01 
01 



Parameters of the HO basis 



Constraints 



TYPE 



TYPE 



TYPE 


TYPE 




0.50 
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MULTCONSTR 
SPINCONSTR 

PRINT- ITER 
EALLMINMAX 

REVIEWFILE 

REVIEW 

RECORDFILE 

REPLAYFILE 

WDDDSAFILE 

CDULOMFILE 

COULOMSAVE 

RESTART 

READ_WDDDS 

NILSSONPAR 

EXECUTE 
ALL_DONE 



1 

-12. 



0.01 42. 1 

0.00 0. 

Output-file parameters 

1 

0. 

Files 



HFODD . REV 
2 

HFODD . REC 
HFODD . REP 
WOODS . WFN 
HFODD . COU 









-1.175 -0.247 



Starting the iteration 



-1.175 -0.352 11.170 11.170 6.280 
Calculate 



Terminate 
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5.2.4. File dyl52-d.dat 



This file (dyl52-d.dat) contains the input data for the code HFODD. 
Only keywords introduced after version (1.60r) are included here. 



General data 



ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

PHASESPACE 

DIASIM_NEU 
DIASIM_PRO 
DIASIG.NEU 
DIASIG_PRO 

COULOMFILE 
COULOMSAVE 



0.0001 



0.01 3 



00 



00 



PARTICLES 

00 00 
PARTICLES 

00 00 
PARTICLES 

00 00 
PARTICLES 

00 00 



00 



00 



Configurations 
HOLES 
00 00 

HOLES 
00 00 
HOLES 
00 00 00 

HOLES 
00 00 00 
— Files 



00 
00 



TYPE 



TYPE 



TYPE 


TYPE 




HFODD . COU 



EXECUTE 
ALL_DONE 



Calculate 



Terminate 
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5.2.5. File dyl52-e.dat 



This file (dyl52-e.dat) contains the input data for the code HFODD. 
Values of all parameters are identical to the default values assumed by 
the code when no input data are provided. 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

PHASESPACE 

SKYRME-SET 

SKYRME-STD 

HBAR20VR2M} 

SPIN.ORBIT 

LANDAU 

LANDAU-SAT 



86 66 
50 
0.0 
1 

0.0 3 




SKM* 

10 

20.73620941 

120.0 120.0 

0.0 0.0 0.0 0.0 0.0 0.0 



Interaction 



-1.0 

EVE_SCA_TS RHD 

1. 1. 
ODD_SCA_TS SPI 

1. 1. 
EVE_SCA_PM RHD 

1. 1. 
ODD_SCA_PM SPI 

1. 1. 

G_SCALING 

1.0 1.0 

INI_FERMI 

-8.0 

INI_DELTA 



1.0 -1.0 

RHOD 
1. 1. 

SPID 
1. 1. 

RHOD 
1. 1. 

SPID 
1. 1. 



8.0 



LPR 
1. 1. 

LPS 
1. 1. 

LPR 
1. 1. 

LPS 
1. 1. 



TAU 
1. 1. 

CUR 
1. 1. 

TAU 
1. 1. 

CUR 
1. 1. 



SCU 
1. 1. 

KIS 
1. 1. 

SCU 
1. 1. 

KIS 
1. 1. 



DIV 
1. 1. 

ROT 
1. 1, 

DIV 
1. 1, 

ROT 
1. 1, 
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FIXDELTA_N 
FIXDELTA_P 
PAIRNFORCE 
PAIRPFORCE 
PAIR_FORCE 
PAIRNINTER 
PAIRPINTER 
PAIR_INTER 
CUTOFF 

SIMPLEXY 

SIGNATUREY 

PARITY 

ROTATION 

TIMEREVERS 

TSIMPLEX_Y 

TSIMPLEXES 

TSIMPLEX3D 

PAIRING 

HFB 

PHNONE_NEU 
PHNONE_PRO 
VACSIM_NEU 
VACSIM_PRO 



1.0 
1.0 
1.0 



1.0 






-200.0 0.0 1.0 

-200.0 0.0 1.0 

-200.0 0.0 1.0 

-200.0 0.16 1.0 

-200.0 0.16 1.0 

-200.0 0.16 1.0 
60.0 



1 
1 

-1 
1 



-1 

1 

1 





1 
1 



PARTICLES 
00 

PARTICLES 
00 

SIMP SIMM 

43 43 
SIMP SIMM 



Symmetries 



Configurations 
HOLES 

00 
HOLES 

00 
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PHSIMP. 
PHSIMP. 
VACPAR. 
VACPAR. 
PHPARI. 
PHPARI. 
VACSIG. 
VACSIG. 
PHSIGN. 
PHSIGN. 
DIANON. 
DIANON. 
DIASIM. 
DIASIM. 
DIAPAR. 
DIAPAR. 
DIASIG. 
DIASIG. 



NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 
NEU 
PRO 



33 33 
PARTICLES 

00 00 
PARTICLES 

00 00 
PARP PARM 

44 42 
PARP PARM 

32 34 
PARTICLES 

00 00 
PARTICLES 
00 00 
PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
16 16 17 17 
PARTICLES 

00 00 
PARTICLES 

00 00 
PARTICLES 
02 

PARTICLES 
02 

PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 

02 02 
PARTICLES 
02 02 



00 



00 



02 



00 



00 



02 



00 



00 



02 



02 



01 



01 



HOLES 
00 00 

HOLES 
00 00 



HOLES 
00 00 

HOLES 
00 00 



HOLES 
00 00 
HOLES 

00 00 
HOLES 

01 
HOLES 

01 
HOLES 

01 01 
HOLES 

01 01 

HOLES 
01 01 

HOLES 
01 01 

HOLES 
01 01 

HOLES 
01 01 



00 
00 



01 
01 



TYPE 


TYPE 


TYPE 



TYPE 



TYPE 



TYPE 



TYPE 


TYPE 




Numerical data 



MAX MULTIP 



MAX_SURFAC 
MAX_MAGNET 
COULOMB 
SLOW-DOWN 






80 
0.5 





79 
0.5 



0.25 
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SLOW-PAIR 
EPS_HERMIT 
OPTI_GAUSS 
GAUSHERMIT 

BASIS_SIZE 
HOMEGAZERO 
SURFAC.PAR 
SURFAC_DEF 
SURFAC_DEF 

OMEGAY 

OMISDY 

OMEGA_XYZ 

OMEGA_RTP 

OMISD_RTP 

OMEGA.TURN 

MULTCONSTR 

SURFCDNSTR 

SPINCDNSTR 

SPICON.XYZ 

SPICON.OME 

PRINT- ITER 
PRINT-MOME 



0.5 



l.OOE-14 



18 18 



32 



Parameters of the HO basis 



15 301 800. 



1.2 



86 



66 1.23 
0.61 
0.10 







Constraints 



0.50 
0.00 

0.00 0.00 0.00 
0.00 0.00 0.00 
0.00 0.00 0.00 



0.01 42. 1 
0.00 00. 



0. 


.00 


0. 





0. 


.00 


0. 





0. 


.00 


0. 





0. 


.00 


0. 





0. 


.00 


0. 






Output-file parameters 

1 
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PRINT-INTR 
EALLMINMAX 
EQUASI_MAX 

REVIEWFILE 

REVIEW 

RECDRDFILE 

RECORDSAVE 

REPLAYFILE 

WDDDSAFILE 

CDULOMFILE 

COULOMSAVE 

RESTART 

CONT.PAIRI 

CONT_OMEGA 

READ_WDDDS 

NILSSONPAR 

EXECUTE 
ALL_DONE 



0. 



1 
1 

-12. 
10.0 

HFODD . REV 
2 

HFODD . REC 
1 

HFODD . REP 
WOODS . WFN 
HFODD . COU 












-1.175 -0.247 



Files 



Starting the iteration 



-1.175 -0.352 11.170 11.170 6.280 
Calculate 



Terminate 
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5.2.6. File tbl51-a.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

VACSIG.NEU 
VACSIG_PRO 
PHSIGN_NEU 
PHSIGN_PRO 
DIASIG.NEU 
DIASIG_PRO 

SLOW-DOWN 

PRINT- ITER 

REVIEWFILE 
COULOMFILE 
COULOMSAVE 

RESTART 



86 65 
2000 
0.0001 
5 

0.01 3 
20 



PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
15 16 17 17 

PARTICLES 
00 00 23 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 00 00 



Conf igurat ions 



HOLES 
00 22 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 00 00 
Numerical data 



0.5 0.5 



Output-file parameters 

1 

Files 



tbl51-a.rev 



tbl51 . cou 



TYPE 



TYPE 




Starting the iteration 



Calculate 



OMEGAY 
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REPLAYFILE 
RECORDFILE 
EXECUTE 

RESTART 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 



0.50 

tbl51-a-0500.rec 
tbl51-a-0500.rec 

1 

0.55 

tbl51-a-0500.rec 
tbl51-a-0550.rec 

0.60 

tbl51-a-0550.rec 
tbl51-a-0600.rec 

0.65 

tbl51-a-0600.rec 
tbl51-a-0650.rec 

0.70 

tbl51-a-0650.rec 
tbl51-a-0700.rec 

0.75 
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REPLAYFILE 
RECDRDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



tbl51-a-0700.rec 
tbl51-a-0750.rec 

0.80 

tbl51-a-0750.rec 
tbl51-a-0800.rec 



ALL_DONE 



Terminate 



5.2. Listings of the test input data files 
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5.2.7. File tbl51-b.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

VACSIG.NEU 
VACSIG_PRO 
PHSIGN_NEU 
PHSIGN_PRO 
DIASIG.NEU 
DIASIG_PRO 

SLOW-DOWN 

PRINT- ITER 

REVIEWFILE 
COULOMFILE 
COULOMSAVE 

RESTART 



86 65 
2000 
0.0001 
5 

0.01 3 
5 



PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
15 16 17 17 

PARTICLES 
00 00 24 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 00 00 



Conf igurat ions 



HOLES 
00 22 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 00 00 
Numerical data 



0.5 0.5 



Output-file parameters 

1 

Files 



tbl51-b . rev 



tbl51 . cou 



TYPE 



TYPE 




Starting the iteration 



Calculate 



OMEGAY 
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REPLAYFILE 
RECORDFILE 
EXECUTE 

RESTART 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 



0.50 

tbl51-b-0500.rec 
tbl51-b-0500.rec 

1 

0.55 

tbl51-b-0500.rec 
tbl51-b-0550.rec 

0.60 

tbl51-b-0550.rec 
tbl51-b-0600.rec 

0.65 

tbl51-b-0600.rec 
tbl51-b-0650.rec 

0.70 

tbl51-b-0650.rec 
tbl51-b-0700.rec 

0.75 



5.2. Listings of the test input data files 



75 



REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



tbl51-b-0700.rec 
tbl51-b-0750.rec 

0.80 

tbl51-b-0750.rec 
tbl51-b-0800.rec 



ALL_DONE 



Terminate 
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5. TEST RUNS 



5.2.8. File tbl51-c.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

VACSIG.NEU 
VACSIG_PRO 
PHSIGN_NEU 
PHSIGN_PRO 
DIASIG.NEU 
DIASIG_PRO 

SLOW-DOWN 

PRINT- ITER 

REVIEWFILE 
COULOMFILE 
COULOMSAVE 

RESTART 



86 65 
2000 
0.0001 
5 

0.01 3 
5 



PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
15 16 17 17 

PARTICLES 
00 00 23 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 24 00 

PARTICLES 
00 00 00 00 



Conf igurat ions 



HOLES 
00 22 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 23 00 

HOLES 
00 00 00 00 
Numerical data 



0.5 0.5 



Output-file parameters 

1 

Files 



tbl51-c . rev 



tbl51 . cou 



TYPE 
0+10 

TYPE 




Starting the iteration 



Calculate 



OMEGAY 



5.2. Listings of the test input data files 
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REPLAYFILE 
RECORDFILE 
EXECUTE 

RESTART 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

SLOW-DOWN 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



0.50 

tbl51-c-0500.rec 
tbl51-c-0500.rec 

1 

0.55 

tbl51-c-0500.rec 
tbl51-c-0550.rec 

0.95 0.95 
0.60 

tbl51-c-0550.rec 
tbl51-c-0600.rec 

0.65 

tbl51-c-0600.rec 
tbl51-c-0650.rec 

0.70 

tbl51-c-0650.rec 
tbl51-c-0700.rec 



78 
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SLOW-DOWN 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



0.5 0.5 
0.75 

tbl51-c-0700.rec 
tbl51-c-0750.rec 

0.80 

tbl51-c-0750.rec 
tbl51-c-0800.rec 



ALL_DONE 



Terminate 
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5.2.9. File tbl51-d.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

VACSIG.NEU 
VACSIG_PRO 
PHSIGN_NEU 
PHSIGN_PRO 
DIASIG.NEU 
DIASIG_PRO 

SLOW-DOWN 

PRINT- ITER 

REVIEWFILE 
COULOMFILE 
COULOMSAVE 

RESTART 



86 65 
2000 
0.0001 
5 

0.01 3 
5 



PPSP PPSM PMSP PMSM 

22 22 21 21 
PPSP PPSM PMSP PMSM 
15 16 17 17 

PARTICLES 
00 00 23 00 

PARTICLES 
00 00 00 00 

PARTICLES 
00 00 24 00 

PARTICLES 
00 00 00 00 



Conf igurat ions 



HOLES 
00 22 00 00 

HOLES 
00 00 00 00 

HOLES 
00 00 23 00 

HOLES 
00 00 00 00 
Numerical data 



0.5 0.5 



Output-file parameters 

1 

Files 



tbl51-d.rev 



tbl51 . cou 



TYPE 
0-10 

TYPE 




Starting the iteration 



Calculate 



OMEGAY 
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REPLAYFILE 
RECORDFILE 
EXECUTE 

RESTART 

SLOW-DOWN 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



0.50 

tbl51-d-0500.rec 
tbl51-d-0500.rec 

1 

0.95 0.95 
0.50 

tbl51-d-0500.rec 
tbl51-d-0500.rec 

0.55 

tbl51-d-0500.rec 
tbl51-d-0550.rec 

0.60 

tbl51-d-0550.rec 
tbl51-d-0600.rec 

0.65 

tbl51-d-0600.rec 
tbl51-d-0650.rec 



5.2. Listings of the test input data files 
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OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 

SLOW-DOWN 

OMEGAY 

REPLAYFILE 

RECORDFILE 

EXECUTE 

OMEGAY 
REPLAYFILE 
RECORDFILE 
EXECUTE 



0.70 

tbl51-d-0650.rec 
tbl51-d-0700.rec 

0.5 0.5 
0.75 

tbl51-d-0700.rec 
tbl51-d-0750.rec 

0.80 

tbl51-d-0750.rec 
tbl51-d-0800.rec 



ALL_DONE 



Terminate 
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5. TEST RUNS 



5.2.10. File lal32-a.dat 



This file (lal32-a.dat) contains the input data for the code HFODD. 

It defines two consecutive runs that lead to a planar solution in 132La. 



NUCLIDE 



ITERATIONS 



ITERAT_EPS 



MAXANTIDSC 



75 57 



5000 



0.000001 



Nuclide 



Convergence 



Skyrme interaction 



SKYRME-STD 



1 



Configurations 

VACSIG.NEU PPSP PPSM PMSP PMSM 

19 19 19 18 
VACSIG_PRO PPSP PPSM PMSP PMSM 

15 15 14 13 

HO basis 



BASIS.SIZE 



SURFAC.PAR 



SURFAC_DEF 



SURFAC_DEF 



MAX_MAGNET 



REVIEWFILE 



REPLAYFILE 



RECORDFILE 



OMEGAY 



MULTCDNSTR 



10 286 800. 



75 57 1.230 



2 0.000 



4 0.000 



3 



lal32-a.rev 



lal32-a.rec 



lal32-a.rec 



0.000 



2 0.5 7.0 



Output data 



File names 



Constraints 



5.2. Listings of the test input data files 
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MULTCONSTR 

RESTART 

EXECUTE 

SIMPLEXY 

SIGNATUREY 

PARITY 

TSIMPLEX3D 

VACPAR_NEU 
VACPAR_PRO 

OMEGA.XYZ 
OMEGA_TURN 
MULTCONSTR 
MULTCONSTR 

RESTART 
CONT_OMEGA 

EXECUTE 
ALL_DONE 



2 2 0.5 7.0 



Starting point 

— Calculate 
Symmetries 





1 





Configurations 



PARP PARM 

38 37 

PARP PARM 

30 27 



Constraints 

0.1443375673 0.1443375673 0.1443375673 1 
1 



2 0.5 7.0 
2 2 0.5 7.0 

1 









Starting point 



Calculate 
Terminate 



84 



5. TEST RUNS 



5.2.11. File lal32-b.dat 



This file (lal32-b.dat) contains the input data for the code HFODD. 

It defines two consecutive runs that lead to a planar solution in 132La. 



NUCLIDE 



ITERATIONS 



ITERAT_EPS 



MAXANTIDSC 



75 57 



5000 



0.000001 



Nuclide 



Convergence 



Skyrme interaction 



SKYRME-STD 



1 



Configurations 

VACSIG.NEU PPSP PPSM PMSP PMSM 

19 19 19 18 
VACSIG_PRO PPSP PPSM PMSP PMSM 

15 15 14 13 

HO basis 



BASIS.SIZE 



SURFAC.PAR 



SURFAC_DEF 



SURFAC_DEF 



MAX_MAGNET 



REVIEWFILE 



REPLAYFILE 



RECORDFILE 



OMEGAY 



MULTCDNSTR 



10 286 800. 



75 57 1.230 



2 0.000 



4 0.000 



3 



lal32-b.rev 



lal32-b.rec 



lal32-b.rec 



0.000 



2 0.5 7.0 



Output data 



File names 



Constraints 



5.2. Listings of the test input data files 
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MULTCONSTR 

RESTART 

EXECUTE 

SIMPLEXY 

SIGNATUREY 

PARITY 

TSIMPLEX3D 

VACPAR_NEU 
VACPAR_PRO 

OMEGA.XYZ 
OMEGA_TURN 
MULTCONSTR 
MULTCONSTR 

RESTART 
CONT_OMEGA 

EXECUTE 
ALL_DONE 



2 2 0.5 7.0 1 

Starting point 



Calculate 

Symmetries 



1 

10 

Configurations 

PARP PARM 

38 37 
PARP PARM 
30 27 

Constraints 

0.0000000 0.1767766953 0.1767766953 1 

1 

2 0.5 7.0 

2 2 0.5 7.0 

Starting point 

1 


Calculate 

Terminate 
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5. TEST RUNS 



5.2.12. File pb208-a.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

PHASESPACE 

SKYRME-SET 
SKYRME-STD 

ROTATION 

VACSIG_NEU 
VACSIG_PRO 

OPTI.GAUSS 
GAUSHERMIT 

BASIS.SIZE 
SURFAC_PAR 
SURFAC_DEF 
SURFAC_DEF 

OMEGAY 
MULTCONSTR 



126 82 
500 

0.0000001 
5 

0.01 3 
5 

64 

SLY4 
1 



64 





54 



42 



42 



Interaction 





Symmetries 



Configurations 

PPSP PPSM PMSP PMSM 

29 29 34 34 
PPSP PPSM PMSP PMSM 
22 22 19 19 

Numerical data 



54 54 



Parameters of the HO basis 

12 455 800. 
126 82 1.23 
2 0.0 

4 0.0 

Constraints 

0.00 
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0. 10 



EALLMINMAX 

REVIEWFILE 

REVIEW 

RECORDFILE 

REPLAYFILE 

COULOMFILE 

COULOMSAVE 

EXECUTE 
BASIS.SIZE 

RESTART 

EXECUTE 

ALL_DONE 



-70. 10. 

pb208-a.rev 
2 

pb208-a.rec 
pb208-a.rec 
pb208-a. cou 
1 1 



0.0 
Output-file parameters 



Files 



26 3654 
1 



800. 



Calculate 



Terminate 
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5. TEST RUNS 



5.2.13. File snl20-a.dat 



General data 



NUCLIDE 

ITERATIONS 

ITERAT_EPS 

MAXANTIOSC 

PING-PONG 

CHAOTIC 

SKYRME-SET 
SKYRME-STD 

INI_FERMI 
INI_DELTA 
PAIR_FORCE 

ROTATION 

PAIRING 

HFB 

VACSIG_NEU 
VACSIG.PRO 

BASIS_SIZE 
SURFAC_PAR 
SURFAC_DEF 
SURFAC_DEF 



70 50 
500 

0.0000001 
5 

0.01 3 
5 

SLY4 

1 



-8.0 -8.0 
1.0 1.0 

-285.88 893.375 1.0 




1 
1 



Interaction 





Pairing 



Symmetries 



Configurations 

PPSP PPSM PMSP PMSM 

22 22 13 13 
PPSP PPSM PMSP PMSM 
12 12 13 13 

Parameters of the HO basis 

14 680 800. 



70 

2 

4 



50 1.23 
0.0 
0.0 



5.2. Listings of the test input data files 
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OMEGAY 
MULTCONSTR 

EALLMINMAX 

REVIEWFILE 

REVIEW 

RECORDFILE 

REPLAYFILE 

COULOMFILE 

COULOMSAVE 

EXECUTE 
ALL DONE 



Constraints 

0.00 

2 0.10 0.0 

Output-file parameters 

-70. 10. 

Files 

snl20-a.rev 
2 

snl20-a.rec 
snl20-a.rec 
snl20-a. cou 
1 1 

Calculate 

Terminate 
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Subject Index 

acronym, 12, 15, 45 
adiabatic, 50, 52, 53 
after variation, 13 
alignment, 26, 53 

angular frequency, 7, 19, 26, 34, 35, 42, 45, 

51-55 

angular momentum, 26, 34, 35, 42, 45, 46, 
53, 55 

average neutron pairing gap, 56 

BCS blocking, 9 
BCS equation, 11, 16 
BCS pairing, 11, 19, 21, 34 
before variation, 13 

canonical state, 11 

Cartesian axes, 46 

Cartesian components, 34 

Cartesian directions, 31, 32, 44, 47, 56 

center-of-mass, 12, 13, 44, 45 

chaotic, 11, 51 

charge block, 11 

charge-parity block, 25 

charge-parity-signature block, 28 

charge-simplex block, 11, 23, 55 

comment, 7, 44, 49 

compatibility, 19, 21, 35 

compiler, 48, 49 

conserved symmetries, 11, 51, 52 
constant potential, 16 
constraint, 30, 34-37, 46, 51, 54 
convergence, 9, 31, 44, 51, 53 
coordinate-space calculations, 56 
Coulomb direct energy, 46 
Coulomb direct potential, 40 
Coulomb exchange energy, 44, 46 
Coulomb field, 30, 40 
Coulomb file, 40, 41, 51 
Coulomb Green function, 30 
Coulomb interaction, 56 
coupling constant, 12, 14, 16, 44, 45 
CPU time, 19, 47-49, 51 
cranking term, 34 



crossing bands, 52 

cutoff energy, 11, 18, 32, 56 

data line, 7, 8, 12, 39, 40 
data set, 7, 39, 46 
data value, 7 

default value, 7-12, 41, 50, 51 
deformation, 44, 45 
deformed harmonic oscillator, 32 
density matrices, 11, 18 
density-dependent part, 15 
density-dependent term, 56 
diabatic blocking, 11, 22, 24, 26, 27, 29, 50, 
51 

diabatic-blocking data, 11 
diagonalization, 11, 55 
divergence, 10, 11, 51, 53-55 
doubly-magic nucleus, 55 

eigenstate, 11, 55 
elementary charge, 44 
end of input data file, 7 
energy, 45 

energy density functional, 34, 36 
equivalent-spectrum, 56 
Euler angles, 45 
execution, 49 

Fermi energy, 16, 56 

first iteration, 16, 42, 45, 55, 56 

form- factor, 17, 56 

(jf-factor, 45 

Gauss-Hermite integration, 31, 47, 56 
Gauss-Hermite quadrature, 44 
generalized spin-orbit, 12, 13 

Hermite polynomials, 31 

HFB calculations, 17, 56 

HFB equation, 21 

HFB pairing, 11, 16, 19, 21, 42 

HO basis, 32, 42-44, 47, 49, 51, 53, 56 

HO quanta, 44, 47 

HO state, 32, 47 



92 



Subject Index 



hole state, 11, 24, 26, 29 

initial mean field, 9, 56 

input item, 7, 50, 51 

input parameter, 7, 41, 47, 50, 51, 55 

intrinsic frame, 38, 45, 46, 55 

intrinsic spin, 53, 56 

item, 7 

iteration, 9-11, 26, 30-32, 35, 38-40, 42, 44, 

45, 51-56 
iteration number, 10 

Kerman-Onishi condition, 55 
keyword line, 7, 8 
kinetic energy, 46 
kinetic-energy operator, 13 

Landau parameters, 13 
linear momentum, 13 

magic numbers, 55 
magnetic moment, 30, 38, 46 
mass, 44 

mass distribution, 54, 55 

mean-field Hamiltonian, 11 

memory, 47, 49 

mixing, 16, 31 

moment of inertia, 46 

multipole moment, 26, 30, 36, 38, 45 

neutron skin, 56 

Nilsson asymptotic state, 45 

Nilsson labels, 46 

Nilsson potential, 9, 32, 43, 51-53, 56 
normalization factor, 36 
nuclear surface, 32, 44, 51 
number of iterations, 9 
number of neutrons, 11 

occupied, 23, 25-28, 46, 53 
occupied state, 22, 23, 25, 55 
odd and odd-odd nuclei, 9 
odd or odd-odd nucleus, 19 
one-body center-of-mass, 13 

pairing, 16-18, 21, 31, 46, 51-53, 55, 56 
pairing energy, 45, 46 
pairing gap, 16, 17, 46, 56 



pairing strength, 16 
pairing window, 11, 56 
parameter set, 12 

parity, 19, 25, 26, 28, 30, 45, 51-55 
parity block, 25, 46, 54, 55 
parity-signature block, 28, 29, 45, 46, 53, 54 
parity-signature configuration, 29, 55 
particle numbers, 9 
particle state, 11, 24, 26, 29 
particle-hole channel, 31 
particle-hole configuration, 11, 45 
particle-hole excitation, 22, 23, 25, 26, 28, 
55 

particle-number symmetry, 56 
particle-particle channel, 56 
Pauli exclusion principle, 22, 23, 25 
phase space, 9 

ping-pong, 10, 11, 51, 53, 54 
planar solution, 53, 55 
principal axes, 32 
principal-axes frame, 45 

quadrupole moment, 26, 45, 51, 54 
quasiparticle energies, 38 
quasiparticle properties, 38 
quasiparticle state, 18, 47, 56 

radius, 45 

rearrangement energy, 46 
recommended value, 11, 50, 51 
record file, 39, 40, 42 
reference configuration, 22, 23, 25, 28 
replay file, 39, 42 
review file, 39 
root-mean-square radii, 56 
Routhian, 27, 38, 45, 46, 52, 53 

saturation density, 17, 56 
self-consistent solution, 10, 55 
signature, 19, 28, 51-53 
simplex, 19, 21, 23, 24, 45, 51-53, 55 
simplex block, 24, 45, 46 
Simpson integration, 30 
single-particle configuration, 22, 52, 54 
single-particle energy, 9, 46 
single-particle observables, 10 
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single-particle properties, 38, 39 

single-particle state, 23, 45, 47 

Skyrme force, 12-14, 44 

Skyrme functional, 9, 44, 45 

solid harmonics, 30 

spatial coordinates, 56 

spherical basis, 55, 56 

spherical code, 56 

spherical ground state, 55, 56 

spherical harmonics, 36 

spin, 26, 35, 45, 46 

spin-orbit, 12, 46 

stability, 9, 10, 46 

stability energy, 45 

standard FORTRAN input, 7 

standard FORTRAN output, 44 

starting point, 56 

superdeformed, 7, 50, 51 

surface multipole moment, 30, 36-38 

tensor, 12 

termination of the HF iteration, 9 
theorem of self-consistent symmetries, 55 
tilted-axis-cranking, 50, 53 
time-even contribution, 45, 46 
time-even coupling constant, 15 
time-even potential, 31 
time-odd contribution, 45, 46 
time-odd coupling constant, 13, 15 
time-odd potential, 9, 31 
time-reversal, 19, 21, 51-53 
traditional spin-orbit, 12 
triaxial, 53-55 
two-body center-of-mass, 13 

vector processor, 47 

wave function, 11, 40 
Woods-Saxon file, 40, 42 

zero-range density-dependent interaction, 56 
zero-range density-dependent pairing, 16, 17 
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Index of variables 

ALPHAR, 32, 33 
ASKEDA, 35 
ASKEDI, 35 

CNILSN, 43 
CNILSP, 43 

DELFIN, 16 
DELFIP, 17 
DELINI, 16, 42 
DNILSN, 43 
DNILSP, 43 

ECUTOF, 11, 18 
EFFSAT, 14 
EMAXAL, 38 
EMAXQU, 38 
EMINAL, 38 
ENECUT, 32 
EPSHER, 31 
EPSITE, 9, 51 
EPSPNG, 10, 11, 51 

FACTGN, 16 
FACTGP, 16 
FCHOMO, 32 
FERINI, 16, 42 
FILCOU, 40, 41 
FILREC, 39 
FILREP, 39, 42 
FILREV, 39 
FILWOO, 40 
FURMAX, 30, 40, 41 

GOXAN, 13 
GOPLAN, 13 
GIXAN, 13 
GIPLAN, 13 

HBANIX, 43 
HBANIY, 43 
HBANIZ, 43 
HBMINP, 13 
HBMSAT, 14 



LCRAY, 48 
ICONTI, 9, 39, 42, 43 
ICOULI, 40, 41, 51 
ICOULO, 40, 41, 51 
IDEFIN, 16, 17 
IDEFIP, 17 
IFLAGA, 35 
IFLAGI, 35 
IFLAGQ, 36 
IFLAGS, 36, 37 
IMOVAX, 35, 42, 45, 55 
IN_FIX, 9, 11 
INNUMB, 32 
INTRIP, 38 
lOCONT, 42, 45 
lOPTGS, 31 

IPAHFB, 9, 11, 16-19, 21, 34, 42 

IPAIRI, 9, 11, 16, 19, 21-23, 25, 28, 34, 35 

IPARTY, 19, 22, 25 

IPCONT, 16, 42 

IPRLN, 38 

IPRLP, 38 

IPRLT, 38 

IPRMID, 38 

IPRSTA, 38 

IPRSTO, 38 

IREAWS, 9, 31, 40, 42, 43 
IREVIE, 39 

IROTAT, 9, 19, 21, 34, 35 

ISIGNY, 19, 21, 23, 25, 28 

ISIMPY, 19, 21-23, 25, 28, 34, 35, 46 

ISIMTX, 19, 21, 34, 35 

ISIMTY, 19, 21, 34, 35 

ISIMTZ, 19, 21, 34, 35 

ISTAND, 12 

ITILAX, 34, 35 

ITIREV, 19, 21 

ITISAX, 34, 35 

IWRIRE, 39, 40 

IZ_FIX, 9, 11 

IZHPEV, 48 

IZNUMB, 32 
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KETA_J, 12, 44 
KETA_M, 12, 13, 44 
KETA_W, 12, 13, 44 
KETACM, 12, 13, 44 
KHFLIG, 29 
KHFLIM, 24, 26, 27 
KHFLIQ, 26, 27 
KHFLIZ, 22 
KHNONE, 22 
KHPMSM, 28, 29 
KHPMSP, 28, 29 
KHPPSM, 28, 29 
KHPPSP, 28, 29 
KHSIMM, 23 
KHSIMP, 23 
KHSIQM, 25 
KHSIQP, 25 
KOFLIG, 29 
KOFLIM, 24 
KOFLIQ, 26, 27 
KOFLIZ, 22 
KPFLIG, 29 
KPFLIM, 24, 26, 27 
KPFLIQ, 26, 27 
KPFLIZ, 22 
KPNONE, 22 
KPPMSM, 28, 29 
KPPMSP, 28, 29 
KPPPSM, 28, 29 
KPPPSP, 28, 29 
KPSIMM, 23 
KPSIMP, 23 
KPSIQM, 25 
KPSIQP, 25 
KVASIG, 28 
KVASIM, 23 
KVASIQ, 25 

LAMBDA, 32, 33, 36, 37 
LANODD, 13, 14 
LDBASE, 32, 47 



MIU, 32, 33, 36, 37 



NDBASE, 47 
NDMAIN, 47 



NDSTAT, 11, 47 
NDXHRM, 47 
NDYHRM, 47 
NDZHRM, 47 
NILDAT, 43 
NLIMIT, 32, 47 
NMACON, 30 
NMAPRI, 30 
NMUCON, 30 
NMUCOU, 30 
NMUPRI, 30 
NOITER, 9, 10, 51 
NOSCIL, 32, 47 
NSICON, 30 
NSIPRI, 30 
NUCHAO, 11, 51 
NULAST, 9, 51 
NUMBSP, 11, 47, 55 
NUMCOU, 30, 40, 41 
NUMETA, 30, 40, 41 
NUPAHO, 22, 23, 25, 28, 29 
NUPING, 10 
NXHERM, 31, 47 
NYHERM, 31, 47 
NZHERM, 31, 47 

OMEGAX, 34, 35, 42 
OMEGAY, 34, 35, 42 
OMEGAZ, 34, 35, 42 
OMEPHI, 34 
OMERAD, 34 
OMETHE, 34 
OMISOX, 34, 35 
OMISOY, 34, 35 
OMISOZ, 34, 35 

POWERN, 17 
POWERP, 17, 18 
POWERT, 17, 18 
PRHO_N, 17 
PRHO_P, 17, 18 
PRHO_T, 17, 18 
PRHODN, 17 
PRHODP, 17 
PRHODT, 17 
PRHOSN, 17 
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PRHOSP, 18 
PRHOST, 18 

QASKED, 36 

ROPARM, 32 
RHOSAT, 14 

SASKED, 36, 37 
SCUR_M, 16 
SCUR_P, 16 
SCUR_S, 15 
SCUR_T, 15 
SDIV_M, 15 
SDIV_P, 15 
SDIV_S, 14, 15 
SDIV_T, 14 
SKIS^, 16 
SKIS_P, 16 
SKIS_S, 15 
SKIS_T, 15 
SKYRME, 12 
SLOWEV, 30, 31 
SLOWOD, 30, 31 
SLOWPA, 31 
SLPR_M, 15 
SLPR_P, 15 
SLPR_S, 14 
SLPR_T, 14 
SLPS_M, 16 
SLPS_P, 16 
SLPS_S, 15 
SLPS_T, 15 
SRHO^'I, 15 
SRHO_P, 15 
SRHO_S, 14 
SRHO_T, 14, 15 
SRHODM, 15 
SRHODP, 15 
SRHODS, 14 
SRHODT, 14 
SROT_M, 16 
SROT_P, 16 
SROT_S, 15 
SROT_T, 15 
SSCU_M, 15 



SSCU_P, 15 
SSCU_S, 14 
SSCU_T, 14 
SSPLM, 16 
SSPLP, 16 
SSPLS, 15 
SSPLT, 15 
SSPIDM, 16 
SSPIDP, 16 
SSPIDS, 15 
SSPIDT, 15 
STAU_M, 15 
STAU_P, 15 
STAU_S, 14 
STAU_T, 14 
STIFFA, 35 
STIFFI, 35 
STIFFQ, 36 
STIFFS, 36, 37 

WOJNP, 12, 13 
WOPINP, 12, 13 

XO_LAN, 13 
X1_LAN, 13 
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ALLJDONE, 7 

BASIS_SIZE, 32, 47 

CHAOTIC, 11 
CONT_OMEGA, 42 
CONT_PAIRI, 42 
COULOMB, 30 
COULOMFILE, 40 
COULOMSAVE, 40 
CUTOFF, 18 

DIANON_NEU, 22 
DIANON_PRO, 22 
DIAPAR_NEU, 22, 26 
DIAPAR_PRO, 27 
DIASIG_NEU, 29 
DIASIG_PRO, 29 
DIASIM_NEU, 24 
DIASIM_PRO, 24 

EALLMINMAX, 38 
EPS_HERMIT, 31 
EQUASLMAX, 38 
EVE_SCA_PM, 15 
EVE_SCA_TS, 14, 45 
EXECUTE, 7, 8, 39, 50, 51 

FIXDELTA_N, 16 
FIXDELTA_P, 17 

G_SCALING, 16 
GAUSHERMIT, 31, 47 

HBAR20VR2M, 13 
HFB, 21 

HOMEGAZERO, 32 

INLDELTA, 16 
INLFERMI, 16 
ITERAT_EPS, 9 
ITERATIONS, 9 

LANDAU, 13 
LANDAU-SAT, 14 



MAX_MAGNET, 30 
MAX_MULTIP, 30 
MAX_SURFAC, 30 
MAXANTIOSC, 9 
MULTCONSTR, 36 

NILSSONPAR, 43 
NUCLIDE, 9 

ODD_SCA_PM, 16 
ODD_SCA_TS, 15, 45 
OMEGA_RTP, 34 
OMEGA_TURN, 35 
OMEGA_XYZ, 34, 42, 44 
OMEGAY, 34 
OMISO_XYZ, 34 
OMISOY, 34 
OPTLGAUSS, 31, 44, 47 

PAIR_FORCE, 17 
PAIRJNTER, 18 
PAIRING, 21 
PAIRNFORCE, 17 
PAIRNINTER, 17, 18 
PAIRPFORCE, 17 
PAIRPINTER, 18 
PARITY, 19 

PHASESPACE, 11, 47, 48 
PHNONE_NEU, 22 
PHNONE_PRO, 22 
PHPARLNEU, 25 
PHPARLPRO, 25 
PHSIGN_NEU, 28 
PHSIGN_PRO, 29 
PHSIMP_NEU, 23 
PHSIMP_PRO, 23 
PING-PONG, 10 
PRINT-INTR, 38 
PRINT-ITER, 38 
PRINT-MOME, 38 

READ_WOODS, 42 
RECORDFILE, 39 
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RECORDSAVE, 40 
REPLAYFILE, 39 
RESTART, 42 
REVIEW, 39 
REVIEWFILE, 39 
ROTATION, 19 

SIGNATUREY, 19 
SIMPLEXY, 19 
SKYRME-SET, 12 
SKYRME-STD, 12, 13, 44 
SLOW-DOWN, 30 
SLOW-PAIR, 16, 31 
SPICON_OME, 35 
SPICON_XYZ, 35 
SPIN_ORBIT, 12, 13 
SPINCONSTR, 35 
SURFAC_DEF, 32, 43 
SURFAC_PAR, 32 
SURFCONSTR, 36 

TIMEREVERS, 19 
TSIMPLEX3D, 21 
TSIMPLEX_Y, 21 
TSIMPLEXES, 21 

VACPAR_NEU, 25 
VACPAR_PRO, 25 
VACSIG_NEU, 28 
VACSIG_PRO, 28 
VACSIM_NEU, 23 
VACSIM_PRO, 23 

WOODSAFILE, 40 
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CGEMM, 48 
CPU_TIME, 49 
CPUTIM, 49 

DCLOCK, 49 

ETIME, 49 

F02AXE, 48 

GETTIM, 49 

PARAMS, 12 
PREDEF, 7 

REVIEW, 39 

SETBAS, 44 

TIME, 49 

ZHPEV, 48 
ZHPEVX, 48 
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dyl52-a.dat, 50, 51, 57 
dyl52-a.out, 50 
dyl52-b.dat, 50, 51, 58 
dyl52-b.out, 50 
dyl52-c.dat, 50, 51, 61 
dyl52-c.out, 50 
dyl52-d.dat, 50, 51, 64 
dyl52-d.out, 50 
dyl52-e.dat, 50, 51, 65 
dyl52-e.out, 50 
dyl52e-dat, 51 

hfodd.f, 47, 49 

lal32-a.dat, 50, 53, 55, 82 
lal32-a.out, 50 
lal32-b.dat, 50, 53, 55, 84 
lal32-b.out, 50 

pbl28-a.dat, 50 
pb208-a.dat, 55, 86 
pb208-a.out, 50 

snl20-a.dat, 50, 56, 88 
snl20-a.out, 50 

tb-151-d.dat, 53 
tbl51-a.dat, 50, 52, 53, 70 
tbl51-a.out, 50 
tbl51-b.dat, 50, 53, 73 
tbl51-b.out, 50 
tbl51-c.dat, 50, 53, 76 
tbl51-c.out, 50 
tbl51-d.dat, 50, 52, 79 
tbl51-d.out, 50 



